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1.  Introduction 


Upon  the  request  of  the  Department  of  State  (DOS)  in  conjunction  with  the  US 
Agency  for  International  Development  (USAID),  the  Department  of  Defense  (DOD) 
may  be  called  to  respond  to  global  disasters  (hurricanes,  earthquakes,  tsunamis, 
flood,  volcanic  eruption,  epidemic,  fire,  etc.)  to  conduct  foreign  humanitarian  assis¬ 
tance  and  disaster  relief  Humanitarian  Assistance/Disaster  Relief  (HA/DR)  opera¬ 
tions  to  relieve  or  reduce  human  suffering,  disease,  hunger,  or  privation.1  HA/DR 
operations  are  often  characterized  by  rapidly  changing  circumstances  and  asymmet¬ 
ric  information,  where  there  is  substantial  pressure  to  quickly  provide  relief  supplies 
and  material  to  an  affected  population  overseas.  If  the  DOD  is  called  upon,  the  US 
Transportation  Command  (USTRANSCOM)  is  typically  responsible  for  planning 
and  executing  HA/DR  aerial  delivery  operations  by  providing  assistance  within  a 
few  days  in  the  form  of  food  and  water  (consumable  aid).1  In  addition  to  the  pro¬ 
vision  of  consumable  aid,  the  DOD  is  often  tasked  to  assist  with  other  humanitar¬ 
ian  services,  such  as  sanitation,  health  care,  nonconsumable  aid  (clothing,  bedding, 
etc.),  emergency  shelter,  and  support  to  critical  infrastructure.1  These  latter  HA/DR 
services  are  not  required  as  quickly  as  consumable  aid. 

Within  the  72  h  immediate  response  authority2  of  a  foreign  humanitarian  disaster, 
there  lacks  an  effective,  safe  ability  to  provide  food  and  water  directly  to  individuals 
or  small  groups  within  the  affected  population  via  ground  delivery  operations.  Nor¬ 
mal  HA/DR  distribution  channels  are  blocked  or  nonexistent,  and  first  responders 
from  non-governmental  organizations  often  cannot  reach  the  affected  population 
because  of  impassible  roads,  inoperable  airports  and  seaports,  and  logistical  bottle¬ 
necks.  Since  disasters  can  render  normal  ground  lines  of  communication  or  trans¬ 
portation  unusable  for  long  periods  of  time,  aerial  delivery  is  the  only  mechanism 
for  providing  immediate  aid.  This  method  safely  delivers  and  disperses  food  and 
water  directly  to  the  affected  population,  thereby  reducing  or  possibly  even  elimi¬ 
nating  the  need  for  on-ground  support  in  the  disaster  area.  Further,  aerial  delivery 
allows  USAID  and  the  affected  state  to  conduct  needs  assessments  and  establish  the 
infrastructure  required  to  provide  sustained  humanitarian  assistance  via  land,  air,  or 
sea.1 

Humanitarian  aid  distribution  to  the  affected  population  by  means  of  aerial  delivery 
has  some  drawbacks.  For  example,  aerial  delivery  does  not  ensure  that  the  most 
vulnerable  people  will  have  access  to  this  airdropped  consumable  aid  (since  nobody 


1 


is  necessarily  coordinating  the  distribution  on  the  ground).  However,  we  emphasize 
that  aerial  delivery  is  the  first  means  to  reach  the  affected  population,  but  it  should 
be  not  maintained  for  long  periods  of  time.  In  other  situations,  humanitarian  aid 
delivery  to  affected  areas  includes  landing  and  giving  it  to  people  in  charge  of  on¬ 
ground  distribution.  In  this  report,  we  do  not  consider  this  latter  situation  but  only 
aerial  delivery  operations  (such  as  precision  airdrop). 

To  support  the  strategic  logistics  planning  and  cost  assessment  of  military  HA/DR 
aerial  delivery  operations  for  the  DOD,  a  decision  support  model  is  needed  to  facil¬ 
itate  both  responsive  and  efficient  logistics  planning  of  the  military  HA/DR  aerial 
delivery  supply  chain  network.  An  optimal  supply  network  design  is  a  function  of 
the  humanitarian  aid  to  be  procured,  transported,  stored,  and  delivered  across  the 
multistage  supply  chain. 

1.1  Problem  Definition 

The  US  Army  Natick  Soldier  Research,  Development  and  Engineering  Center 
(NSRDEC),  with  support  from  USTRANSCOM  and  many  other  partners  across 
the  DOD  and  DOS,  has  initiated  development  of  the  next  generation  HA/DR  aerial 
delivery  capability  to  enable  airdrops  to  occur  directly  over  the  affected  popula¬ 
tion.  The  direct  aerial  delivery  of  HA/DR  consumable  aid  (food  and  water)  through 
airdrop  operations  avoids  the  need  for  ground  supply  chain  lines,  which  are  often 
damaged  as  a  result  of  the  disaster,  while  providing  HA/DR  aid  within  days  follow¬ 
ing  a  disaster.3 

In  conjunction  with  the  development  of  the  next-generation  HA/DR  airdrop  capa¬ 
bility,  NSRDEC  is  further  developing  the  food  items,  water  items,  and  airdrop  con¬ 
tainers,  and  they  are  investigating  HA/DR  aid  distribution  effectiveness  and  spread 
modeling.  In  addition  to  these  efforts,  NSRDEC  seeks  a  strategic  logistics  plan¬ 
ning  model  for  military  HA/DR  aerial  delivery  operations  in  order  to  understand 
better  the  trade-offs  of  supply  chain  efficiency  and  responsiveness  in  the  event  of  a 
disaster.3 

NSRDEC  would  like  to  identify  the  optimal  number  and  location  of  pre-positioned 
facilities  (distribution  centers)  for  storing  HA/DR  consumable  aid  supplies.  This 
entails  analyzing  current  US  Air  Force  bases  with  consideration  of  the  location 
of  cargo  aircraft  capable  of  responding,  proximity  to  high  disaster  risk  areas,  and 
storage  space  available.  Further,  the  amount  of  consumable  humanitarian  aid  (in- 
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ventory)  to  be  both  procured  and  stored  must  be  optimally  determined  to  meet  the 
estimated  affected  population  size  (consumer  demand).  Moreover,  NSRDEC  de¬ 
sires  a  cost  assessment  capability  (i.e.,  an  analysis  of  potential  cost  savings)  of  the 
entire  military  HA/DR  supply  network  from  the  procurement  of  food  and  water 
items  to  the  actual  HA/DR  aerial  delivery  mission  execution.3 

1.2  Motivation  and  Purpose 

In  this  project,  we  develop  a  Multiple  Criteria  Decision  Analysis  (MCDA)  frame¬ 
work  to  optimize  the  military  HA/DR  aerial  delivery  supply  chain  network.  The 
model  uses  stochastic,  mixed-integer,  weighted  goal  programming  to  optimize  net¬ 
work  design,  logistics  costs,  staging  locations,  procurement  amounts,  and  inven¬ 
tory  levels.  As  a  result,  the  optimization  model  provides  strategic  decision  support 
to  NSRDEC  and  other  HA/DR  decision  makers  for  the  aerial  delivery  of  small 
HA/DR  airdrop  systems:  airdrop  of  small  food  and  water  (consumable  aid)  pack¬ 
ages  directly  onto  affected  disaster  relief  populations. 

The  MCDA  framework  enables  decision  makers  to  explore  the  trade-offs  between 
military  HA/DR  aerial  delivery  supply  chain  efficiency  and  responsiveness  while 
optimizing  across  a  wide  range  of  real-world  probabilistic  scenarios  to  account  for 
the  inherent  uncertainty  in  the  location  and  severity  of  global  humanitarian  disas¬ 
ters  as  well  as  the  amount  of  demand  to  be  met.  Moreover,  the  MCDA  framework 
provides  strategic  decision  support  to  help  answer  the  following  questions: 

•  From  a  set  of  candidate  staging  facilities  around  the  world,  how  many  and 
which  ones  should  be  opened  to  pre-position  and  distribute  consumable  aid 
to  disaster  areas  via  several  aerial  delivery  modes? 

•  How  much  consumable  aid  should  be  procured  from  the  suppliers  as  well  as 
stored  in  inventory  at  both  depots  and  staging  facilities? 

•  How  many  trips  are  necessary  (between  the  different  supply  chain  stages)  for 
each  of  the  different  pre -positioned  aerial  delivery  modes  to  achieve  a  best 
compromise  solution,  trading  off  supply  chain  efficiency  and  responsiveness? 
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2.  Literature  Review 


Supply  chain  management  for  HA/DR  delivery  operations  is  the  process  of  plan¬ 
ning,  implementing,  and  controlling  the  efficient,  cost-effective  flow  and  storage 
of  goods  and  materials,  as  well  as  related  information,  from  the  point  of  origin  to 
the  point  of  consumption  for  the  purpose  of  alleviating  the  suffering  of  vulnerable 
people.4  Logistics  in  the  HA/DR  sector  encompasses  several  traditional  activities, 
such  as  the  procurement,  transportation,  and  warehousing  of  goods  and  services, 
as  well  as  other  specific  activities,  such  as  disaster  preparedness  and  planning.4^6 
Network  design  emphasizes  the  optimal  design  of  the  HA/DR  supply  chain  using 
mathematical  models  and  methods  to  determine  optimal  strategies  and  policies  for 
managing  the  supply  chain.7  Important  performance  criteria  include  supply  chain 
efficiency,  responsiveness,  and  risk.7 

Characteristics  that  bring  additional  complexity  and  unique  challenges  to  HA/DR 
delivery  supply  chain  design  and  management  include  unpredictability  of  demand 
in  terms  of  timing,  geographic  location,  type  of  commodity,  and  quantity  of  com¬ 
modity;  suddenly  occurring  demand  in  very  large  amounts  and  short  lead  times  for  a 
wide  variety  of  supplies;  high  importance  associated  with  response  time;  and  lack  of 
initial  resources  (supply,  people,  technology,  transportation  capacity,  and  money).8 
The  flow  of  resources  coincides  with  4  main  phases  of  disaster  relief:  assessment 
(minimal  resources  are  required  to  identify  what  is  needed),  deployment  (resource 
requirements  ramp  up  to  meet  the  needs),  sustainment  (operations  are  sustained  for 
a  period  of  time),  and  reconfiguration  (operations  are  reduced,  then  terminated).8 

Given  that  HA/DR  delivery  supply  chains  usually  operate  in  highly  uncertain  en¬ 
vironments,  they  must  be  engineered  and  executed  in  shorter  periods  of  time  as 
to  provide  relief  to  the  affected  population  as  soon  as  possible.9  Further,  inventory 
management  in  HA/DR  delivery  supply  chains  is  affected  by  unreliable,  incom¬ 
plete  or  nonexistent  information  about  lead  times,  demand  levels,  and  locations.10 
In  terms  of  distribution  network  configuration,  the  number  and  location  of  distribu¬ 
tion  centers  is  uncertain.  This  makes  cost  assessment  difficult  in  terms  of  planning 
financial  flows. 

For  the  HA/DR  supply  chain  network,  there  are  3  dominating  costs:  supply  costs, 
distribution  costs,  and  inventory  holding  costs.  Unpredictable  demand  patterns  in¬ 
crease  the  complexity  of  relief  organization-supplier  relationships,  making  them 
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more  difficult  to  foster  than  in  the  relatively  stable  demand  environment  of  the 
commercial  supply  chain. 1 1  Further,  supply  procurement  options  generally  cannot 
be  evaluated  before  a  disaster  occurs.  Thus,  it  may  be  difficult  to  control  the  cost 
of  supplies.  Distribution  costs  stem  from  the  need  to  transport  massive  amounts  of 
materials  in  a  very  short  amount  of  time.  Varied  disaster  locations  lead  to  varied 
modes  of  transportation.  The  types  of  inventory  costs  include  inventory  investment, 
inventory  obsolescence  (and  spoilage),  order/setup  costs,  and  holding  (carrying) 
costs.  Inventory  control  for  supply  warehouses  in  the  relief  chain  is  challenging 
because  of  the  higher  variations  in  lead  times,  demands,  and  demand  locations.11 
Coordinating  activities  between  various  agencies  for  HA/DR  is  also  a  challenging 
issue. 12-14 

Relatively  few  studies  have  been  conducted  to  optimize  HA/DR  delivery  opera¬ 
tions  in  terms  of  distribution  center  pre-positioning  and  inventory  staging.  A  sur¬ 
vey  of  decision  aid  models  for  humanitarian  logistics  is  provided  in  Ortuno  et  al.15 
Akkihal16  solved  an  array  of  mixed-integer  linear  program  formulations  to  examine 
the  strategic  impact  of  inventory  pre-positioning  on  delivery  lead  time  of  HA/DR 
operations.  The  model  determines  optimal  locations  for  warehousing  nonconsum¬ 
able  inventories  required  for  initial  deployment  of  aid.  The  objective  of  the  model 
is  to  minimize  the  average  global  distance  from  the  nearest  warehouse  to  a  fore¬ 
casted  homeless  person.  Demand  patterns,  along  with  correlated  variables,  such  as 
population  and  hazard  frequency,  offer  views  of  regional  vulnerability  to  natural 
disasters.16  Similarly,  Balcik  and  Beamon8  considered  facility  location  decisions 
for  a  humanitarian  relief  chain  responding  to  the  quick-onset  disasters,  where  they 
develop  an  optimization  model  that  determines  the  number  and  locations  of  distri¬ 
bution  centers  in  a  relief  network  and  the  amount  of  relief  supplies  to  be  stocked 
at  each  distribution  center  to  meet  the  needs  of  the  affected  population.  The  model 
integrates  facility  location  and  inventory  decisions  while  considering  multiple  item 
types  and  capturing  budgetary  constraints  and  capacity  restrictions.8 

Duran  et  al.17  developed  a  mixed-integer  programming  inventory  location  model  to 
evaluate  the  effect  that  pre-positioning  relief  items  would  have  on  average  relief- 
aid  emergency  response  time.  They  found  the  optimal  number  and  location  of  pre¬ 
positioning  warehouses  given  that  demand  for  relief  supplies  can  be  met  from  both 
pre-positioned  warehouses  and  suppliers.  They  allowed  multiple  HA/DR  events  to 
occur  within  a  replenishment  period,  thus  capturing  the  adverse  effect  of  ware¬ 
house  replenishment  lead  time.  They  also  allowed  the  probability  of  need  for  each 
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item  to  depend  on  both  local  conditions  and  natural  hazard  type.17  Salmeron  and 
Apte18  developed  a  2-stage  stochastic  optimization  model  with  the  goal  of  mini¬ 
mizing  casualties  by  determining  location  and  expansion  decisions  of  assets,  such 
as  warehouses  and  shelters,  in  the  first  stage  and  then  determining  the  best  sup¬ 
porting  logistics  decisions  in  the  second  stage.  The  second-stage  decisions  classify 
the  population  into  those  that  need  emergency  evacuation  (critical),  those  that  need 
commodities  (stay-back),  and  those  that  are  displaced  (transfer).  The  stochastic  op¬ 
timization  formulation  helps  model  the  inherent  uncertainty.18 

Mogilevsky19  developed  the  Disaster  Relief  Airlift  Planner  (DRAP),  which  is  an 
optimization-based  decision- support  tool  that  determines  optimal  routes  to  deliver 
material  given  certain  data  such  as  disaster  location  and  available  airports,  aircraft, 
and  supply  stockpiles.  DRAP  is  formulated  to  minimize  disaster  material  short¬ 
ages  while  preferring  to  choose  routes  that  reduce  transportation  costs  (and  delivery 
times)  based  on  decision-maker  constraints  and  priorities.  The  model  is  also  useful 
for  helping  determine  the  optimal  aircraft  allocation  and  positioning  for  HA/DR  op¬ 
erations.  DRAP  can  be  used  by  logistics  planners  and  decision  makers  to  conduct 
trade-off  analysis  among  routes  with  respect  to  transportation  costs  and  demand 
shortages  in  very  short  time  horizon  logistics  planning.19 

Even  fewer  studies  have  used  multiple  criteria  optimization  approaches  for  improv¬ 
ing  HA/DR  delivery  operations  and  logistics  planning.  As  an  exception,  Park20  de¬ 
veloped  a  multiobjective  decision-making  model  to  incorporate  the  decision  maker(s) 
value  trade-offs  in  the  disaster  relief  resource  allocation  problem.  The  decision  win¬ 
dow  for  resource  allocation  is  the  critical  first  72  h  after  the  initial  damage  assess¬ 
ment  has  been  made.  Value-focused  thinking  was  used  to  capture  the  value  trade¬ 
offs,  and  the  resulting  value  hierarchy  is  optimized  via  a  mathematical  programming 
model  to  solve  the  multiobjective  resource  allocation  problem.20  In  another  study, 
Vitoriano  et  al.21  proposed  a  goal  programming  model  to  provide  decision  support 
to  solve  the  multicriteria  humanitarian  logistics  aid  distribution  problem,  attempt¬ 
ing  to  minimize  costs  and  time  of  response  while  maximizing  equity  of  distribution 
or  reliability  and  security  of  the  operation  routes.21 

In  addition,  Bozorgi-Amiri  et  al.22  developed  a  multiobjective  robust  stochastic  pro¬ 
gramming  approach  for  HA/DR  logistics  under  uncertain  demand,  supply,  and  costs 
(procurement  and  transportation).  Further,  the  model  considers  uncertainty  for  the 
locations  where  those  demands  might  arise  and  the  possibility  that  some  of  the  pre- 
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positioned  supplies  in  the  relief  distribution  center  or  supplier  might  be  partially 
destroyed  by  the  disaster.  The  multiobjective  model  attempts  to  minimize  the  sum 
of  the  expected  value  and  the  variance  of  the  total  cost  of  the  relief  chain  while 
penalizing  the  solution’s  infeasibility  due  to  parameter  uncertainty.  The  model  also 
aims  to  maximize  the  affected  areas’  satisfaction  levels  through  minimizing  the 
sum  of  the  maximum  shortages  in  the  affected  areas.  To  solve  this  bicriteria  prob¬ 
lem,  they  formulated  a  compromise  programming  model  and  solve  it  to  obtain  a 
Pareto-efficient  (nondominated)  compromise  solution.  The  purpose  of  the  model  is 
to  provide  decision  support  on  both  facility  location  and  resource  allocation  in  cases 
of  HA/DR  efforts.22 

Despite  these  recent  studies,  a  critical  gap  in  the  literature  remains  with  assessing 
the  trade-offs  between  supply  chain  efficiency  (i.e.,  total  logistics  costs)  and  supply 
chain  responsiveness  (i.e.,  supply  delivery  time,  demand  fulfillment)  for  aerial  de¬ 
livery  operations  in  the  entire  military  HA/DR  supply  chain  network.  None  of  these 
previous  studies  considered  the  delivery  of  HA/DR  consumable  aid  via  aerial  deliv¬ 
ery  mechanisms.  Further,  none  of  these  previous  studies  considered  the  trade-offs 
of  response  time,  total  cost,  and  amount  of  demand  satisfied.  Therefore,  we  seek  to 
fill  the  gap  in  the  literature  while  providing  military  HA/DR  decision  makers  with 
strategic  decision  support  using  our  MCDA  framework. 

3.  Materials  and  Methods 

In  this  section,  we  first  briefly  describe  the  military  HA/DR  aerial  delivery  process 
flow  for  consumable  aid  distribution.  We  then  briefly  depict  the  military  HA/DR 
aerial  delivery  supply  chain  network.  Last,  we  discuss  the  MCDA  framework  and 
formulate  the  stochastic,  mixed-integer,  weighted  goal  programming  model  used 
for  strategic  military  HA/DR  logistics  planning. 
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3.1  Military  HA/DR  Aerial  Delivery  Logistics  Network 

In  the  event  of  a  natural  disaster,  military  HA/DR  aerial  delivery  planners  and  policy 
makers  must  make  a  series  of  quick  decisions  in  an  effort  to  provide  humanitarian 
aid  to  the  affected  population  as  soon  as  possible.  In  order  to  facilitate  this  tacti¬ 
cal  and  operational-level  decision  making,  strategic  planners  must  fully  understand 
the  process  flow  of  consumable  aid  for  military  HA/DR  aerial  delivery  operations. 
Figure  1  depicts  the  military  HA/DR  aerial  delivery  process  flow. 


Fig.  1  Military  HA/DR  aerial  delivery  process  flow 


In  this  aerial  delivery  process  flow  depicted  in  Fig.  1,  HA/DR  consumable  aid 
is  purchased  from  the  vendor(s)  and  transported  to  (typically  stateside)  HA/DR 
supply  depot(s).  Next,  the  consumable  aid  is  transported  from  the  depot(s)  to  pre¬ 
positioned  HA/DR  staging  locations  around  the  world,  where  the  consumable  aid 
is  stored.  In  the  event  that  a  natural  disaster  occurs  requiring  HA/DR  aerial  deliv¬ 
ery  operations,  the  aerial  delivery  capability  is  immediately  prepared  for  mission 
execution,  and  the  consumable  aid  is  transported  and  air  dropped  to  the  affected 
population. 

In  addition  to  the  process  flow,  Fig.  2  depicts  the  military  HA/DR  aerial  delivery 
supply  chain  network. 
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Supplier  Warehouse  Distribution  Beneficiary 


Fig.  2  Military  HA/DR  aerial  delivery  supply  chain  network 


The  military  HA/DR  aerial  delivery  supply  chain  network  in  Fig.  2  matches  the  pro¬ 
cess  flow  illustrated  in  Fig.  1.  Upon  examining  the  4  stages  of  the  supply  chain  net¬ 
work,  we  see  that  the  vendors  represent  the  suppliers,  the  depots  represent  the  ware¬ 
houses,  the  pre-positioned  storage  facilities  represent  the  distribution  centers,  and 
the  disaster  areas  represent  the  beneficiaries.  We  use  this  multistage  supply  chain 
network  representation  in  Fig.  2  to  develop  the  stochastic,  mixed-integer,  weighted 
goal  programming  optimization  model.  This  MCDA  framework  is  used  as  a  strate¬ 
gic  planning  tool  for  optimizing  the  military  HA/DR  aerial  delivery  supply  chain 
network  design. 

3.2  Multiple  Criteria  Decision  Analysis  Framework 

In  our  MCDA  framework,  we  incorporate  several  modeling  methodologies  to  pro¬ 
vide  a  robust  decision  support  tool  (see  Appendix  A).  First,  we  use  goal  program¬ 
ming  to  allow  decision  makers  to  systematically  explore  and  examine  different 
optimization  problem  criteria.  The  decision  maker  defines  goals  for  each  of  the 
criteria  considered  and  then  evaluates  the  effects  each  of  these  criteria  have  on 
the  global  optimal  solution.  This  methodology  is  particularly  useful  for  strategic 
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planning  when  incorporated  with  goal  priority  weights  determined  by  the  decision 
maker(s).23  Second,  we  use  Design  of  Experiments  (DOE)  to  estimate  the  impact 
of  the  underlying  factors  causing  uncertainty  within  the  system.  We  use  a  23  full- 
factorial  DOE  to  model  the  interactive  effects  between  experimental  design  factors; 
this  approach  helps  identify  robust  alternatives  over  the  set  of  probabilistic  scenar¬ 
ios.23  Third,  we  use  stochastic  optimization  to  incorporate  random  elements  into  the 
model  objective  function  and  data  parameters;  this  method  provides  more  robust  so¬ 
lutions  to  aid  decision  makers  when  optimizing  under  uncertainty.23  In  this  model, 
we  optimize  the  expected  value  of  the  objective  function  (due  to  the  probabilistic 
scenarios)  and  calculate  many  of  the  data  parameters  stochastically. 

Our  MCDA  framework  uses  stochastic,  mixed-integer,  weighted  goal  program¬ 
ming.  We  next  describe  the  optimization  model  assumptions,  sets,  parameters,  de¬ 
cision  variables,  objective,  and  constraints. 

3.2.1  Optimization  Model  Assumptions 

In  this  optimization  model,  we  assume  that  each  delivery  mode  transports  its  maxi¬ 
mum  allowable  capacity  each  trip  and  that  all  consumable  aid  on  board  is  delivered. 
We  assume  that  each  aerial  delivery  mode  is  capable  of  traveling  the  full  distance 
required  each  trip  in  terms  of  flight  crew  and  fuel  (i.e.,  fuel  is  available  en  route,  if 
necessary).  We  also  assume  that  there  is  no  limit  on  the  number  of  hours  flown  for 
each  aerial  delivery  mode.  Finally,  we  assume  that  all  pre-positioned  aerial  delivery 
assets  are  available  to  be  used  (if  necessary,  but  number  of  trips  may  vary)  to  deliver 
consumable  aid  to  the  disaster  area  (aircraft  will  fly  in  parallel,  as  opposed  to  each 
aircraft  waiting  for  the  prior  vehicle  to  return  before  departing  with  aid). 

3.2.2  Optimization  Model  Sets 

The  following  list  represents  the  optimization  model  sets: 

V  -  set  of  consumable  aid  vendors  (suppliers)  with  v  £  V 
D  =  set  of  consumable  aid  depots  (warehouses)  with  d  £  D 
M  =  set  of  candidate  storage  facilities  (distribution  centers)  with  j  £  M 
N  =  set  of  high-risk  disaster  areas  (beneficiaries)  with  i  £  N 
K  =  set  of  aerial  delivery  modes  with  k  £  K 
G  =  set  of  decision-maker  goals  with  g  £  G 
S  =  set  of  disaster  planning  scenarios  with  s  £  S 
B  =  set  of  model  iterations  with  b  £  B 
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3.2.3  Optimization  Model  Parameters 

The  following  list  represents  the  optimization  model  parameters: 

psb  =  probability  of  occurrence  of  disaster  planning  scenario  s  in  iteration  b 
wg  =  decision-maker  weight  for  goal  g 

Disrupts b  =  supply  chain  disruption  effect  on  supply  chain  costs  in  scenario  s  and 
iteration  b 

Foodsb  =  number  of  days  worth  of  food  to  provide  in  scenario  s  and  iteration  b 

Impsb  =  disaster  impact  above  or  below  estimated  demand  in  scenario  s  and  iter¬ 
ation  b 

Derriisb  =  expected  consumable  aid  demand  (pounds)  per  day  at  disaster  area  i  in 
scenario  s  and  iteration  b 

dsij  =  geodesic  distance  (nautical  miles)  from  storage  facility  j  to  disaster  area  i 

ldj  =  geodesic  distance  (nautical  miles)  from  depot  d  to  storage  facility  j 

Spdk  =  average  speed  of  aerial  delivery  mode  k 

LoadDelayk  =  load  time  for  aid  (per  trip)  of  aerial  delivery  mode  k 

tijk  =  travel  time  (per  trip)  from  facility  j  to  disaster  area  i  via  delivery  mode  k 

rdj  =  travel  time  (per  trip)  from  depot  d  to  storage  facility  j 

OMFk  =  operation  and  maintenance  and  fuel  cost  (dollars  per  hour)  of  aerial 
delivery  mode  k 

T Dcostdj  =  transportation  cost  (per  trip)  from  depot  d  to  storage  facility  j 

TScostijk  =  transportation  cost  (per  trip)  from  facility  j  to  disaster  area  i  via  de¬ 
livery  mode  k 

Atd  =  amount  of  aid  (pounds)  transported  (per  trip)  from  depot  d  to  storage  facility 

j 

Atsk  =  amount  of  aid  (pounds)  transported  (per  trip)  from  facility  j  to  area  i  via 
mode  k 
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IHDcostd  =  inventory  holding  cost  (per  pound)  at  depot  d 
IHScostj  =  inventory  holding  cost  (per  pound)  at  storage  facility  j 
Pcostv  -  procurement  cost  (per  pound)  of  consumable  aid  from  vendor  v 
Capj  =  inventory  holding  capacity  (pounds)  at  storage  facility  j 
Cpd  =  inventory  holding  capacity  (pounds)  at  depot  d 
Fcostj  =  fixed  cost  of  opening  storage  facility  j 

Ntddj  =  maximum  number  of  trips  with  consumable  aid  from  depot  d  to  storage 
facility  j 

NtSijk  =  maximum  number  of  trips  with  consumable  aid  from  storage  facility  j  to 
disaster  area  i  via  aerial  delivery  mode  k 

ADMjk  =  number  of  aerial  delivery  modes  k  pre-positioned  at  storage  facility  j 

TGlsb  =  target  goal  for  the  total  aerial  delivery  response  time  in  scenario  s  and 
iteration  b 

TG 2sb  =  target  goal  for  for  total  supply  chain  cost  in  scenario  s  and  iteration  b 
TG 3sb  =  target  goal  for  for  unmet  demand  in  scenario  s  and  iteration  b 

3.2.4  Optimization  Model  Decision  and  Goal  Deviation  Variables 

The  following  list  represents  the  optimization  model  decision  and  goal  deviation 
variables: 

Xj  =  1  if  open  storage  facility  j,  or  0  otherwise 

Yij  =  1  if  disaster  area  i  is  served  by  storage  facility  j,  or  0  otherwise 

Zl]ksb  =  number  of  trips  with  consumable  aid  from  storage  facility  j  to  disaster  area 

i  via  aerial  delivery  mode  k  in  scenario  s  and  iteration  b 

Hdjsb  =  number  of  trips  with  consumable  aid  from  depot  d  to  storage  facility  j  in 
scenario  s  and  iteration  b 

Qdv  =  amount  of  aid  (pounds)  purchased  from  vendor  v  for  storage  at  depot  d 
Invj  =  amount  of  aid  (pounds)  to  store  in  inventory  at  storage  facility  j 
Irid  =  amount  of  aid  (pounds)  to  store  in  inventory  at  depot  d 
P0SgSb  =  positive  deviation  for  goal  g  in  scenario  s  and  iteration  b 
neggsb  =  negative  deviation  for  goal  g  in  scenario  s  and  iteration  b 
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3.2.5  Optimization  Model  Formulation 

The  following  is  a  description  of  the  optimization  model  objective  function  and 
constraints: 


min  EE  Psb{Wg=\pOSg=\tSb  +  W g=2POS g=2,sb  T  W  g=‘.iPO^  g='i,sb )  ,  (1) 

b  s 

subject  to  EEE  tijk^ijksb  EE  f djHdjsb  POSg=\,sb  TGI  ,  (2) 

i  j  k  d  j 

Vs  e  S,b  e  B 

Disrupts b  (EE  T  Dcostdj  Hdjsb 

d  j 

+EEEE  T  S  co  stij  k  Z ij  ksb 

i  j  k  v 

+  ( FcostjXj  +  IHScostjlnVj)  ,  (3) 

j 

+  51  IHDcostdInd  +  55  55  -P  COStyQ  dv^ 

d  v  d 

—  posg=2,sb  =  TG2  Vs  G  S',  b  G  5 

55  DemisbImpsbF oodsb  -  55  55  55  AtSkZijksbADMjk  ,  (4) 

i  i  j  k 

—  posg=3tsb  =  TG3  Vs  e  S',  6  e  -B 


J5  V './  -  1  ViGiV,  (5) 

j 

55  E,  =  X,  Vj  e  M  ,  (6) 

i 

EE  AtskZijksbADMjk  <  Invj  Vj  e  M,  s  e  S,b  e  B  ,  (7) 

i  fc 

Invj  <  CapjXj  Vj  G  M  ,  (8) 

/nnj  >  Xj  Vj  G  M  ,  (9) 

Atd 55  E  Ind  Vd  £  D,  s  <E  S,b  <E  B  ,  (10) 

j 

Ind<Cpd  Vd  e  D  ,  (11) 

5 ~^Qdv<Ind  VdeD,  (12) 

V 
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Hdjsb  <  NtddjXj  Vd  e  D,j  e  M,s  e  S,b  e  B  ,  (13) 

Zijksb  <  NtSijkYij  Vi  £  N,j  £  M,k  £  K,s  £  S,b  £  B  t  (14) 

Atd  E  Hdjsb  ^  EE  AtSkZijksbADMjk  ,  (15) 

d  i  k 

Vj  e  M,s  e  S,be  B 

^’eV,seS,beB  ,  (16) 


d  d  j 

Xj  e  (o,  1}  Vj  e  M 
Yij  e  (0, 1}  VieNjeM 

^ijksb  >o,  integer  Vi  e  N,  j  G  M,  k  &  K,s  £  S,b  £  B 

Hdjsb  >  0,  integer  VdeD}jeM,s£S,beB 

Invj  >  0  Vj  £  M  ,  (17) 

Irid  >0  Vd  £  D 

Qdv  E  0  Vc?  £  d,  v  £  V 

P0SgSb  >0  Vg  £  G,  s  £  S,b  £  B 

neggsb  >0  Vg  £  G,  s  £  S,b  £  B 

The  objective  function  in  Eq.  1  seeks  to  minimize  the  sum  of  the  three  expected 
weighted  goal  deviations  for  target  response  time,  target  budget,  and  target  de¬ 
mand  met,  across  all  probabilistic  disaster  planning  scenarios  s  and  iterations  b.  The 
weights  allow  the  planner  to  prioritize  the  targets.  Goal  constraints  in  Eq.  2  ensure 
target  TGI  is  met  for  the  total  aerial  delivery  response  time  (in  trip  hours)  for  each 
scenario  s  and  iteration  b.  The  amount  by  the  target  is  not  met,  the  posg=liSi ,  >  0, 
is  minimized  in  Eq.  1;  TGI  is  set  to  0  since  the  decision  maker  wishes  to  minimize 
total  aerial  delivery  time,  which  is  a  function  of  replenishment  time  between  depots 
and  storage  facilities  as  well  as  the  final  disaster  response  time  between  storage 
facilities  and  affected  areas. 

Goal  constraints  in  Eq.  3  ensure  the  target  for  total  supply  chain  cost  (TG2)  with 
supply  chain  disruption  factor  is  met  for  each  scenario  s  and  iteration  6;  target  de¬ 
viation  is  captured  in  posg=2,sb ■  The  amount  by  which  the  target  is  not  met,  the 
posg=2,sb  E  0,  is  minimized  in  Eq.  1;  TG2  is  set  to  0  since  the  decision  maker 
wishes  to  minimize  total  supply  chain  cost.  Goal  constraints  in  Eq.  4  ensure  that  the 
amount  of  delivered  consumable  aid  shortage  to  the  affected  population  including 
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disaster  impact  and  food  factors  for  each  scenario  s  and  iteration  b  does  not  exceed 
the  target  TG 3;  the  amount  of  shortage  above  the  target  is  captured  in  posg=z,sb- 
Again,  the  amount  by  which  the  target  is  not  met,  the  posg=s,sb  >  0,  is  minimized 
in  Eq.  1;  TG3  is  set  to  0  since  the  decision  maker  desires  for  the  estimated  demand 
to  be  met. 

The  constraints  in  Eq.  5  ensure  that  each  disaster  area  i  is  served  by  exactly  one 
storage  facility  j.  The  constraints  in  Eq.  6  ensure  that  if  disaster  area  i  is  served 
by  storage  facility  j,  then  storage  facility  j  must  be  opened  for  pre-positioning 
inventory  and  assets  and  that  each  storage  facility  j  can  only  serve  one  disaster  area 
i;  Eqs.  5  and  6  are  set  partitioning  constraints.  From  Eq.  7,  the  amount  of  inventory 
to  store  at  storage  facility  j  must  be  greater  than  or  equal  to  the  sum  of  the  product 
of  the  number  of  trips,  the  amount  of  consumable  aid  transported  per  trip,  and  the 
number  of  assets  pre-positioned  for  every  scenario  s  and  iteration  b.  From  Eq.  8, 
the  inventory  is  held  only  at  opened  storage  facilities,  and  the  amount  of  inventory 
kept  at  the  storage  facility  j  must  not  exceed  its  capacity.  The  constraints  in  Eq.  9 
enforce  that  if  a  storage  facility  j  is  opened,  then  inventory  must  be  stored  there. 

The  constraints  in  Eq.  10  ensure  that  the  amount  of  inventory  to  store  at  each  depot 
d  must  be  greater  than  or  equal  to  the  sum  of  the  product  of  the  number  of  trips  (of 
a  single  C-17  aircraft)  and  the  amount  of  consumable  aid  transported  per  trip  for 
every  scenario  s  and  iteration  b.  In  Eq.  11,  all  depots  are  assumed  to  be  open,  but 
the  amount  of  inventory  kept  at  depot  d  must  not  exceed  its  capacity.  In  Eq.  12,  if 
consumable  aid  is  purchased  from  the  vendors  for  storage  at  depot  d,  then  there  is 
inventory  for  storage  at  depot  d.  The  constraints  in  Eq.  13  ensure  that  the  number 
of  trips  with  consumable  aid  from  each  depot  d  to  each  storage  facility  j  must  not 
exceed  its  maximum;  this  maximum  is  a  function  of  a  target  total  flight  time  budget 
value  set  by  the  decision  maker. 

In  Eq.  14,  the  number  of  trips  with  consumable  aid  from  each  storage  facility  j  to 
each  disaster  area  i  via  each  delivery  mode  k  must  not  exceed  its  maximum;  again, 
this  maximum  is  a  function  of  a  target  total  flight  time  budget  value  by  the  decision 
maker.  The  constraints  in  Eq.  15  ensure  that  the  total  amount  of  consumable  aid 
shipped  from  the  depots  d  to  each  storage  facility  j  must  appropriately  backfill  the 
total  amount  of  aid  delivered  to  the  disaster  areas  i  for  each  scenario  s  and  iteration 
b.  From  Eq.  16,  the  amount  procured  from  each  of  the  vendors  v  for  storage  at  all 
of  the  depots  d  must  be  greater  than  or  equal  to  the  total  amount  of  aid  shipped  out 
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of  depot  d  for  each  scenario  s  and  iteration  b.  In  Eq.  17,  we  have  the  binary,  integer, 
and  non-negativity  decision  variable  constraints. 

3.2.6  Optimization  Model  Scenarios  and  Iterations 

As  part  of  the  multiple  criteria  decision  analysis  framework,  we  incorporated  a  23 
full-factorial  experimental  design  as  a  mechanism  to  help  mitigate  uncertainty  asso¬ 
ciated  with  consumable  aid  demand  as  well  as  the  effects  of  supply  chain  disruption 
on  total  supply  chain  costs.  By  exploring  3  different  design  factors  (each  at  2  levels), 
we  generate  8  different  disaster  planning  scenarios  that  are  combinatorially  applied 
to  each  of  the  disaster  areas  in  the  military  HA/DR  supply  chain  network.  Each  of 
these  disaster  areas  represent  a  unique  model  iteration,  making  sb  scenario-iteration 
pairs.  These  scenarios  and  iterations  are  discussed  further  in  Section  4.2. 

4.  Computational  Experiment  and  Results 

In  this  section,  we  describe  the  computational  experiments  performed  for  NSRDEC 
using  the  MCDA  framework  for  strategic-level  military  HA/DR  supply  chain  net¬ 
work  design  optimization.  The  results  of  this  computational  experiment  provided 
the  decision  maker  with  the  necessary  decision  support  to  assess  the  supply  chain 
efficiency  and  responsiveness  trade-offs  for  military  HA/DR  aerial  delivery  opera¬ 
tions. 

4.1  Data  Collection  and  Analysis 

The  military  HA/DR  aerial  delivery  supply  chain  network  consists  of  4  separate 
stages.  For  the  first  stage,  we  consider  v  =  3  consumable  aid  vendors  (suppliers), 
where  the  supplier  selection  subproblem  consists  of  determining  how  much  aid 
to  purchase  from  each  vendor  v  for  storage  at  each  depot  d.  This  supply  network 
design  uses  a  multiple  sourcing  strategy.7  Given  the  uncertainty  associated  with 
purchase  prices,  we  model  the  procurement  costs  for  each  vendor  v  ( PCostv )  by 
randomly  sampling  from  uniform  distributions  with  minimum  and  maximum  values 
determined  by  the  decision  maker  based  on  subject-matter  expertise. 

For  the  second  stage,  we  consider  d  =  2  consumable  aid  depots  (warehouses)  in 
the  continental  United  States  as  a  means  of  risk  pooling  (risk  pooling  is  the  practice 
of  consolidating  facilities  into  fewer  locations  by  using  a  consolidated  distribution 
strategy);  risk  pooling  can  help  reduce  supply  chain  costs  as  well  as  demand  risk.7 
For  this  stage  in  the  supply  chain  network,  the  2  depots  considered  were  US  Marine 
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Corps  Logistics  Bases  (MCLBs)  in  Albany  (Georgia,  USA)  and  Barstow  (Califor¬ 
nia,  USA).  For  this  distribution  and  inventory  planning  subproblem,  the  model  de¬ 
termines  the  amount  of  aid  to  store  in  inventory  at  depot  d  as  well  as  the  number  of 
trips  with  consumable  aid  from  depot  d  to  storage  facility  j.  We  assume  that  there  is 
one  aerial  delivery  mode  (C-17  aircraft)  available  at  each  depot  for  transportation, 
and  that  the  amount  of  aid  transported  per  trip  from  depot  d  to  storage  facility  j 
(. Atd )  in  scenario  s  and  iteration  b  equals  its  maximum  loading  capacity.  The  in¬ 
ventory  holding  costs  at  each  depot  ( IHDcostd )  and  inventory  holding  capacity  at 
each  depot  (Cpa)  are  modeled  by  randomly  sampling  from  uniform  distributions 
with  minimum  and  maximum  values  again  determined  using  subject-matter  exper¬ 
tise. 

For  the  third  stage  of  the  supply  chain  network,  we  consider  m  —  10  candidate 
storage  facilities  (distribution  centers)  for  pre-positioning  k  =  3  different  aerial 
delivery  assets  (C-17  Globemaster,  C-130  Hercules,  CH-47  Chinook)  and  consum¬ 
able  aid  inventory.  The  10  storage  facilities  considered  were  US  Air  Force  Bases  lo¬ 
cated  in  both  the  United  States  (North  Carolina,  Washington,  Florida,  and  Hawaii) 
as  well  as  overseas  in  foreign  countries  (Germany,  Japan,  Turkey,  South  Korea, 
Kuwait,  and  Honduras).  This  stage  of  the  supply  chain  network  design  uses  a  de¬ 
consolidated  distribution  strategy  given  that  there  are  more  facilities,  each  serving 
a  smaller  demand  region;  this  supply  chain  strategy  can  help  increase  supply  chain 
responsiveness,  reduce  the  risk  of  supply  chain  disruption,  and  reduce  outbound 
transportation  costs.7  For  this  facility  location,  inventory,  and  distribution  planning 
subproblem,  the  model  determines  the  candidate  storage  facilities  j  to  open,  the  dis¬ 
aster  areas  i  served  by  each  facility  j,  the  number  of  trips  with  consumable  aid  from 
facility  j  to  disaster  area  i  using  aerial  delivery  mode  k,  and  the  amount  of  aid  to 
store  in  inventory  at  each  storage  facility  j.  Again,  we  assume  that  the  amount  of  aid 
transported  per  trip  from  facility  j  to  disaster  area  i  via  each  aerial  delivery  mode  k 
(. AtSk )  equals  the  maximum  loading  capacity  of  each  asset.  The  inventory  holding 
costs  at  each  facility  ( IHScostj ),  the  inventory  holding  capacity  at  each  facility 
0 Capj ),  and  the  fixed  cost  of  opening  each  storage  facility  (Fcostj)  are  modeled  by 
randomly  sampling  from  uniform  distributions  with  minimum  and  maximum  val¬ 
ues  again  determined  by  the  decision  maker  using  subject-matter  expertise.  Also, 
the  number  of  aerial  delivery  modes  k  available  for  pre-positioning  at  each  storage 
facility  j  ( ADMjk )  is  determined  by  the  decision  maker. 

To  calculate  the  geodesic  distance  (in  nautical  miles)  between  each  depot  d  and 
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storage  facility  j  ( l dj),  we  used  the  Haversine  formula.24  This  method  estimates 
the  great-circle  distance  between  2  points  on  a  sphere  using  their  longitudes  and 
latitudes.  Equation  18  shows  the  distance  calculation: 


Idj  —2rc  arcsin 


+  cos(0i)  cos (02)  sin2 


Vd  e  DJ  e  M 


(18) 


where  r  is  the  radius  of  the  sphere  (6,371  km  for  Earth),  c  is  the  kilometers  to 
nautical  miles  conversion  factor  (0.539),  0i  and  Ai  are  the  latitude  and  longitude, 
respectively,  of  the  first  point,  and  02  and  A2  are  the  latitude  and  longitude,  respec¬ 
tively,  of  the  second  point. 

In  order  to  mitigate  uncertainty  in  the  transportation  times  and  costs,  we  conducted  a 
Monte  Carlo  simulation  of  1,000  trials  to  compute  the  average  cruise  speed  (knots) 
of  each  aerial  delivery  mode  k  ( Spdk)  and  the  load  delay  time  (hours)  of  each 
aerial  delivery  mode  k  ( LoadDelai/k ).  In  each  trial  of  the  simulation,  both  the  cruise 
speeds  and  load  delay  times  were  randomly  sampled  from  uniform  distributions 
using  minimum  and  maximum  values  based  on  subject-matter  expertise.  Spdk  and 
LoadDelayk  were  computed  as  the  average  values  over  all  trials.  The  travel  time 
(per  one-way  trip)  from  depot  d  to  storage  facility  j  (rdj)  was  calculated  using  Eq. 
19: 

f dj  =  t v~r - ^  LoadDelayk=Cn  Vd  G  DJ  G  M  .  (19) 

bpdk=cu 

As  a  component  of  the  transportation  costs,  there  is  an  operation,  maintenance,  and 
fuel  cost  (in  dollars  per  hour)  for  each  aerial  delivery  asset  k  ( OFMJ ,  which  is 
computed  by  randomly  sampling  from  uniform  distributions  using  minimum  and 
maximum  values  determined  by  the  decision  maker.  Therefore,  the  transportation 
cost  (per  trip)  from  depot  d  to  storage  facility  j  (T  DcostJ  is  calculated  using  Eq. 
20: 


TDcostdj  =  (/  +  OMFk=Cn )  *  rdj  Vd  G  DJ  e  M  ,  (20) 

where  /  is  the  fixed  aircrew  cost  per  hour;  we  use  /  =  250  as  determined  by  the 
decision  maker.  Also,  we  computed  the  maximum  number  of  trips  with  consumable 
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aid  from  depot  d  to  storage  facility  j  ( Ntddj )  using  Eq.  21: 


Ntddj  = 


T 


2  *  r  dj 


Vd,  E  D,j  e  M  , 


(21) 


where  T  is  the  target  total  flight  time  budget  value  between  the  depots  and  facilities 
set  by  the  decision  maker  (an  upper  bound  on  total  flight  hours  available  from  depot 
d).  This  parameter  is  used  to  evaluate  supply  chain  efficiency  and  responsiveness 
trade-offs  in  the  model  sensitivity  analysis.  For  instance,  T  =  120  infers  a  total 
flight  time  budget  of  120  h  (or  5  days). 


In  the  fourth  stage  of  the  supply  chain  network,  we  consider  n  —  8  foreign  dis¬ 
aster  areas  (beneficiaries)  in  separate  high-risk  geographic  regions  using  the  last 
20  years  of  global  disaster  data  from  the  Emergency  Events  Database  (EM-DAT) 
maintained  by  the  Centre  for  Research  on  Epidemiology  of  Disasters  at  the  Uni¬ 
versity  catholique  de  Louvain  and  supported  by  several  organizations  around  the 
world,  to  include  USAID;  EM-DAT  contains  data  on  the  occurrence  and  effects  of 
over  18,000  mass  disasters  in  the  world  from  1900  to  present.25  To  determine  the 
location  of  these  foreign  disaster  areas  representing  the  centroids  of  different  high- 
risk  geographic  regions,  we  performed  both  k-means  clustering  and  hierarchical 
clustering  methods  using  several  variables,  to  include  disaster  frequency,  mortality, 
total  number  affected,  and  total  economic  damages. 

Using  complete  linkage  hierarchical  clustering,  we  determined  8  clusters  was  an  ap¬ 
propriate  number  via  a  given  cut  of  the  associated  dendrogram.  Fork-means  cluster¬ 
ing,  we  identified  the  appropriate  number  of  clusters  via  a  plot  of  the  Within-Group 
Sum  of  Squares  (WGSS)  as  well  as  using  the  results  of  26  different  criteria  (from 
the  NbClust  function  in  R).  The  left  plot  of  Fig.  3  displays  the  WGSS  versus  num¬ 
ber  of  clusters,  where  we  see  the  curve  levels  off  at  8  clusters.  The  right  plot  of 
Fig.  3  shows  a  histogram  of  the  number  of  clusters  chosen  by  each  of  the  criteria, 
which  indicates  a  bimodal  distribution.  Although  the  highest  frequency  of  criteria 
was  3  clusters,  8  clusters  was  also  significant,  as  well  as  more  pragmatic,  for  the 
computational  experiment. 
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Fig.  3  Determining  appropriate  number  of  cluster  for  k-means  clustering 


Both  clustering  methods  identified  the  following  8  foreign  disaster  areas  to  serve 
as  epicenters  of  various  geographic  regions:  Haiti  (Caribbean),  Indonesia  (South- 
Eastern  Asia),  Mexico  (Central  America),  Tanzania  (Eastern  Africa),  India  (South¬ 
ern  Asia),  China  (Eastern  Asia),  Australia  (Australia  and  New  Zealand),  and  Peru 
(South  America).  Figure  426  illustrates  the  locations  of  the  2  consumable  aid  depots 
(squares),  10  candidate  storage  facilities  (diamonds),  and  8  disaster  areas  (stars). 


Fig.  4  Location  of  depots,  storage  facilities,  and  disaster  areas 
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Further,  we  used  the  EM-DAT  disaster  data  with  a  sample  size  of  n  =  2,  797  disas¬ 
ters  (from  1994  to  2014)  to  fit  an  ordinary  least  squares  multiple  regression  model 
to  estimate  the  total  number  of  affected  people  per  year.  To  satisfy  the  assumption 
of  normally  distributed  residuals  (error  terms),  we  performed  a  Box-Cox  transfor¬ 
mation  on  the  response  variable,  which  suggested  a  natural  log  power  transforma¬ 
tion  (A  =  0).  We  then  performed  5-fold  cross-validation  (with  best  subset  selection 
within  each  of  the  training  sets)  to  select  the  best  predictors  for  the  model;  this  re¬ 
sulted  in  the  following  5  variables:  year,  number  of  occurrences  (occur),  number 
of  deaths  (death),  total  monetary  damages  in  thousands  of  dollars  (damage),  and 
population  (pop).  There  was  no  multicollinearity  present  in  the  model  (variance  in¬ 
flation  factors  <  5),  but  the  model  diagnostics  suggested  nonconstant  variance  of 
the  residuals;  the  Modified  Levene’s  Test,  Breusch  Pagan  Test,  and  Special  White’s 
Test  all  confirmed  heteroscedasticity.  Therefore,  we  calculated  robust  standard  er¬ 
rors  to  compensate  for  the  unknown  pattern  of  nonconstant  error  variance,  giving 
more  accurate  p-values  for  the  predictors.  Also,  there  was  no  temporal  or  spatial 
autocorrelation  in  the  data  used  in  the  regression  model.  The  final  regression  model 
results  (K2  =  0.21,  p-value  <  0.0001)  are  displayed  in  Table  1. 


Table  1  Regression  model  results  for  estimated  affected  population  per  year 


Parameter 

Estimate 

Robust  Standard  Error 

t-value 

p-value 

(Intercept) 

-90.316 

29.362753 

-3.08 

0.002 

year 

0.04854 

0.0146549 

3.31 

0.001 

occur 

0.29683 

0.028824 

10.3 

0.000 

death 

1.7E-05 

8.331E-06 

2.06 

0.040 

damage 

1.3E-08 

9.36E-09 

1.41 

0.158 

pop 

-3E-09 

1.13E-09 

-2.68 

0.007 

From  this  regression  model,  we  used  Eq.  22  to  compute  the  expected  consumable 
aid  demand  per  day  at  each  high-risk  disaster  area  i  per  each  scenario  s  and  iteration 
b  ( Derriisb ).  For  each  disaster  planning  scenario  s,  we  assume  in  this  computational 
experiment  that  demand  emanates  from  only  one  high-risk  disaster  area  that  repre¬ 
sents  the  epicenter  of  the  geographic  region  (i.e.,  there  is  not  simultaneous  demand 
at  multiple  disasters  areas).  Thus,  each  disaster  area  i  maps  to  each  iteration  b, 
and  we  ensure  robustness  in  the  supply  chain  network  design  by  optimizing  over 
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sb  =  (8)  (8)  =  64  scenario-iteration  combinations. 


Demisb 


0XT~)  f  ~F  ~F g ^ \ 

: — - — — - —(FoodConsume  +  Water  Consume)  , 

AvgOccuri 

Vi  e  N, s  e  S,b  e  B 


(22) 


where  Af  fecti  is  the  log  of  the  total  number  of  affected  per  year  for  each  disas¬ 
ter  area  i  estimated  using  the  regression  model,  AvgOccuri  is  the  average  (over 
last  20  years)  of  the  total  number  of  disaster  occurrences  per  year  for  each  disas¬ 
ter  area  i,  FoodConsume  is  the  expected  amount  (pounds)  of  food  consumed  per 
affected  person  per  day  (we  used  0.8  lb  as  determined  by  the  decision  maker),  and 
W aterC onsume  is  the  expected  amount  of  water  (pounds)  consumed  per  affected 
person  per  day  (we  used  1.84  lb  as  determined  by  the  decision  maker). 

Using  the  EM-DAT  data,  Table  2  depicts  the  HA/DR  consumable  aid  demand  cal¬ 
culations,  where  Year  is  the  current  year,  Occur  represents  AvgOccuri,  Death  is 
the  average  (over  last  20  years)  of  the  total  number  of  deaths  per  year  caused  by 
disasters,  Damage  is  the  average  (over  last  20  years)  of  the  total  monetary  damages 
in  thousands  of  dollars  per  year  caused  by  disasters,  Pop  is  the  current  population 
for  each  of  the  disaster  areas,  exp  {Affect}  is  the  average  total  number  of  affected 
people  per  year,  and  Demand  is  the  estimated  expected  consumable  aid  demand 
(pounds  per  day)  per  disaster  occurrence  for  each  of  the  disaster  areas  representing 
the  different  worldwide  high-risk  geographic  regions  ( Demisb ). 


Table  2  Estimated  expected  consumable  aid  demand 


Disaster  Area 

Year 

Occur 

Death 

Damage 

Pop 

exp  {Affect} 

Demand 

Australia 

2014 

6 

57 

1,426,790 

23,558,900 

7,983 

9,231 

China 

2014 

60 

7,107 

18,010,076 

1,365,590,000 

1,398,028,739 

36,418,118 

Haiti 

2014 

5 

11.951 

429,331 

10,413,211 

6,690 

10,243 

India 

2014 

36 

5,641 

1,819,772 

1,246,670,000 

1,556,967 

116,250 

Indonesia 

2014 

20 

9,197 

1,273,830 

252,164,800 

306,658 

69,943 

Mexico 

2014 

11 

273 

1,731,776 

119,713,203 

25,377 

14,300 

Peru 

2014 

9 

340 

48,764 

30,814,175 

16,427 

1,2413 

Tanzania 

2014 

6 

411 

0 

44,928,923 

7,028 

8,334 
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To  calculate  the  geodesic  distance  (in  nautical  miles)  between  each  storage  facility  j 
and  disaster  area  i  ( dsij ),  we  used  the  Haversine  formula  again.  Equation  23  shows 
the  distance  calculation: 


Vi  G  N,j  G  M 


The  travel  time  (per  one-way  trip)  from  storage  facility  j  to  disaster  area  i  via  aerial 
delivery  mode  k  ( Ujk )  was  calculated  using  Eq.  24: 

ds  •  • 

t^k  =  +  LoadDelayk  Vi  G  N,j  G  M,  k  G  K  .  (24) 

Spdk 

Further,  the  transportation  cost  (per  trip)  from  storage  facility  j  to  disaster  area  i  via 
aerial  delivery  mode  k  ( TScostijk )  is  calculated  using  Eq.  25: 

TScostijk  =  (/  +  OMFk )  *  tijk  Vi  G  N,  j  G  M,  k  G  K  .  (25) 

In  addition,  we  computed  the  maximum  number  of  trips  with  consumable  aid  from 
storage  facility  j  to  disaster  area  i  via  aerial  delivery  mode  k  ( NtSijk )  using  Eq.  26: 

R 

Nts^k  =  - — - —  Vi  G  N,j  G  M,  k  G  K  ,  (26) 

4  *  tijk 

where  R  is  the  target  total  flight  time  budget  value  between  the  facilities  and  the 
disaster  areas  set  by  the  decision  maker.  Again,  this  parameter  is  used  to  evalu¬ 
ate  supply  chain  efficiency  and  responsiveness  trade-offs  in  the  model  sensitivity 
analysis.  For  instance,  R—  72  infers  a  total  flight  time  budget  of  72  h  (or  3  days). 

4.2  Disaster  Planning  Scenarios 

As  previously  noted,  we  incorporated  a  23  full-factorial  DOE  as  a  mechanism  to 
develop  disaster  planning  scenarios  to  better  account  for  the  uncertainty  within  the 
military  HA/DR  aerial  delivery  supply  chain  network.  The  3  experimental  design 
factors  included  Impsb,  Foodsb,  and  Disrupt sb.  Impsb  represents  a  disaster  impact 
factor  that  adjusts  the  estimated  expected  consumable  aid  demand  in  each  scenario 
s  and  iteration  b  up  or  down  by  25%  to  account  for  the  deviation  from  expected 
impact.  Foodsb  represents  a  days-of-food  factor  that  adjusts  the  estimated  expected 
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consumable  aid  demand  in  scenario  s  and  iteration  b  by  the  number  of  day’s  worth 
of  food  to  provide  for  each  affected  person;  this  factor  is  also  varied  at  2  levels 
(7  or  14  days  worth  of  food).  These  first  2  DOE  factors  address  uncertainty  in¬ 
volving  the  amount  of  beneficiary  demand,  affecting  both  supply  chain  efficiency 
and  responsiveness.  Disrupts b  represents  a  supply  chain  disruption  impact  factor 
to  better  account  for  uncertainty  in  supply  chain  costs  for  each  scenario  s  and  itera¬ 
tion  b,  while  also  increasing  the  robustness  of  the  supply  chain  network  design.  We 
assume  that  supply  chain  disruptions  (such  as  aerial  delivery  weather  delays  and 
maintenance  issues)  will  affect  supply  chain  processes  associated  with  transporta¬ 
tion  and  storage  of  consumable  aid,  which  is  reflected  in  the  overall  supply  chain 
costs.  Therefore,  this  factor  increases  the  total  supply  chain  costs  by  0%  (low  risk) 
and  11%  (high  risk)  (the  average  effect  of  disruptions  is  an  1 1%  growth  in  cost27). 

From  the  8  disaster  planning  scenarios  coupled  with  the  8  model  iterations  (one  it¬ 
eration  per  each  disaster  area),  the  model  optimizes  over  64  total  scenario-iteration 
combinations  to  determine  the  best  compromise  solution  that  balances  supply  chain 
efficiency  and  responsiveness.  We  calculated  the  probability  of  occurrence  of  sce¬ 
nario  s  in  iteration  b  (psb)  using  Eq.  27: 

psh  =  -DisPropb  Vs  e  S,b  G  B  ,  (27) 

s 

where  DisPropb  represents  the  proportion  of  historical  (over  the  last  20  years) 
disaster  frequency  for  each  iteration  b,  which  equals  the  quotient  of  the  average 
number  of  disaster  occurrences  per  year  for  each  disaster  area  (mapping  to  a  unique 
iteration)  divided  by  the  sum  of  the  average  number  of  disaster  occurrences  per  year 
over  all  disaster  areas. 

In  Table  3,  we  highlight  the  23  full-factorial  experimental  design  factors  (and  levels) 
used  to  generate  the  8  disaster  planning  scenarios  for  the  computational  experiment. 

Again,  these  8  disaster  planning  scenarios  shown  in  Table  3  are  applied  to  each  of 
the  8  model  iterations,  representing  a  total  of  64  different  scenario-iteration  combi¬ 
nations  employed  in  the  optimization  model  as  a  means  to  help  mitigate  uncertainty 
in  the  system.  This  method  helps  provide  a  robust  best-compromise  solution  to  the 
stochastic,  mixed-integer  weighted  goal  programming  problem. 
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Table  3  The  23  full-factorial  DOE  factors  for  disaster  planning  scenarios 


Scenario 

Disaster  Impact 

Days  of  Food 

Supply  Chain  Disruption 

1 

0.75 

7 

1 

2 

1.25 

7 

1 

3 

0.75 

14 

1 

4 

1.25 

14 

1 

5 

0.75 

7 

1.11 

6 

1.25 

7 

1.11 

7 

0.75 

14 

1.11 

8 

1.25 

14 

1.11 

4.3  Model  Implementation  and  Solutions 

To  implement  the  MCDA  framework  for  solving  the  multicriteria  military  HA/DR 
aerial  delivery  supply  chain  network  design  optimization  problem  under  uncer¬ 
tainty  and  provide  decision  support  to  NSRDEC,  we  leveraged  the  General  Al¬ 
gebraic  Modeling  System  (GAMS),28  Microsoft  Excel,  and  Microsoft  Visual  Basic 
for  Applications  (VBA)  platforms  (see  Appendix  B).  In  particular,  we  used  GAMS 
v.23.9.3  with  IBM  ILOG  CPLEX  12.4.0.1  to  solve  the  stochastic,  mixed-integer 
weighted  goal  programming  model,  and  we  used  Excel/VBA  to  create  an  auto¬ 
matic,  user-friendly  interface  with  the  decision  maker  for  model  input  and  analysis 
of  model  output.  The  following  best  compromise  solution  to  the  computational  ex¬ 
periment  was  solved  in  0.499  s  on  a  Lenovo  Thinkpad  W510  laptop  with  an  Intel  i7 
CPEi  (1.6  GHz)  and  8.00  GB  of  random-access  memory. 

In  the  following  solution,  the  decision  maker  heavily  weighted  the  third  goal  to 
place  significant  priority  on  meeting  the  estimated  expected  consumable  aid  de¬ 
mand  across  all  high-risk  disaster  areas.  The  solution  depends  upon  the  pseudo¬ 
random  number  generator  seed  in  GAMS.  Table  4  highlights  which  storage  facili¬ 
ties  were  opened,  the  amount  of  pre-positioned  inventory  to  store,  and  the  number 
of  required  trips  by  each  aerial  delivery  model. 

Given  that  each  storage  facility  must  serve  exactly  one  disaster  area,  the  results  in 
Table  4  indicate  that  8  of  the  10  candidate  storage  facilities  are  opened  for  pre¬ 
positioning  both  aerial  assets  and  consumable  aid  inventory.  The  US  Air  Force 
Bases  in  Kuwait  and  Honduras  should  not  be  opened  as  part  the  optimal  network. 
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Table  4  Results  for  storage  facility  locations,  inventory,  and  trips 


Storage  Facility  Name 

Opened 

Inventory 

Number  of  Trips 

Ramstein  Air  Base  (Germany) 

1 

665,600 

4 

2 

0 

Charlotte  ANG  Base  (North  Carolina,  USA) 

1 

832,000 

3 

5 

0 

McChord  Air  Force  Base  (Washington,  USA) 

1 

896,000 

4 

5 

0 

Hickam  Air  Force  Base  (Hawaii,  USA) 

1 

769,600 

3 

2 

1 

Yokota  Air  Base  (Japan) 

1 

640,000 

3 

5 

0 

Incirlik  Air  Base  (Turkey) 

1 

588,800 

3 

4 

0 

Osan  Air  Base  (South  Korea) 

1 

563,200 

2 

4 

0 

Ali  A1  Salem  Air  Base  (Kuwait) 

0 

0 

0 

0 

0 

Soto  Cano  Air  Base  (Honduras) 

0 

0 

0 

0 

0 

MacDill  Air  Force  Base  (Florida,  USA) 

1 

435,200 

4 

C-17a 

7 

C-130a 

0 

CH-47a 

“C-17  Globemaster,  C-130  Hercules,  CH-47  Chinook. 


In  addition,  Table  5  highlights  the  amount  of  inventory  to  store  at  each  of  the  state¬ 
side  MCLBs  as  well  as  the  breakdown  of  the  amount  procured  from  each  of  the  3 
consumable  aid  vendors.  These  results  in  Table  5  indicate  that  both  depots  in  the 
supply  chain  network  are  used  to  store  inventory  stateside,  but  it  is  clear  that  the 
depot  in  Albany,  Georgia  (USA),  is  used  more  heavily.  Further,  only  one  of  the  3 
vendors  supplies  consumable  aid  to  both  depots.  Figure  5  shows  the  number  of  trips 
via  C-17  from  each  of  the  stateside  depots  to  the  opened  storage  facilities;  this  is 
the  average  over  all  scenario-iteration  combinations.  As  anticipated,  there  are  no 
trips  from  either  depots  to  the  US  Air  Forces  Bases  in  Kuwait  (Ali  A1  Salem)  and 
Honduras  (Soto  Cano)  because  those  facilities  were  not  opened. 


Table  5  Results  for  depot  inventory  and  vendor  procurement 


Depot  Name 

Inventory 

Amount  Procured  from  Vendors 

MCLB  Albany  (Georgia,  USA) 
MCLB  Barstow  (California,  USA) 

14,208,000 

2,496,000 

3,072,000  5,568,000  5,568,000 

2,496,000  0  0 

Vl°  V2a  V3a 

“VI  =  vendor  1,  V2  =  vendor  2,  V3  =  vendor  3. 
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Number  of  Trips  from  Depots  to  Storage  Facilities 
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Storage  Facility 
Ramstein  Air  Base 
Charlotte  Air  National  Guard  Base 
McChoid  Air  Force  Base 
Hidcam  Air  Force  Base 
Yokota  Air  Base 
Indrlik  Air  Base 
Osan  Air  Base 
Ali  Al  Salem  Air  Base 
Soto  Cano  Air  Base 
MacDill  Air  Foroe  Base 


Marine  Corps  Logistics  Base  Albany  Marine  Corps  Logistics  Base  Barstow 


Fig.  5  Results  for  number  of  trips  from  depots  to  storage  facilities 


Figure  5  also  indicates  that  the  distribution  of  trips  from  the  depot  in  Barstow,  Cal¬ 
ifornia  (USA),  to  all  of  the  opened  storage  facilities  are  relatively  uniform.  Further, 
the  greatest  number  of  trips  are  from  the  depot  in  Albany,  Georgia  (USA),  to  the 
storage  facility  at  Charlotte  Air  National  Guard  (ANG)  Base  in  Charlotte,  North 
Carolina  (USA);  this  makes  sense  because  these  2  locations  are  relatively  close  in 
terms  of  geodesic  distance. 

In  addition,  Table  6  depicts  which  global  disaster  areas  (representing  high-risk  ge¬ 
ographic  regions)  are  served  by  which  opened  storage  facilities. 

From  these  best  compromise  solution  results,  Fig.  626  provides  a  visual  illustration 
of  the  optimal  military  HA/DR  aerial  delivery  supply  chain  network  design. 

Upon  investigating  the  trade-offs  of  supply  chain  efficiency  versus  responsiveness 
for  this  best-compromise  solution  to  the  multiple-criteria  problem,  the  goal  devia¬ 
tion  variables  provided  critical  information  about  response  time,  cost,  and  demand 
unmet  for  every  scenario  s  and  iteration  b.  In  particular,  this  supply  chain  network 
design  provided  an  average  response  time  (across  all  disaster  areas)  of  roughly  6 
days,  an  average  total  supply  chain  cost  (across  all  stages  of  the  supply  chain  net¬ 
work)  of  roughly  $153  billion,  and  an  average  total  demand  shortage  (across  all 
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disaster  areas)  of  roughly  47  million  pounds.  Moreover,  this  supply  chain  network 
design  provided  identical  median  values  for  the  response  time  and  cost,  but  the  me¬ 
dian  total  demand  shortage  equaled  0;  in  fact,  there  were  only  8  of  the  64  scenario- 
iteration  combinations  (13%)  where  the  demand  was  unmet. 


Table  6  Results  for  disaster  area  and  storage  facility  assignments 


Disaster  Area  (Geographic  Region)  Servicing  Storage  Facility 


Haiti  (Caribbean) 

Indonesia  (Southeast  Asia) 
Mexico  (Central  America) 
Tanzania  (Eastern  Africa) 

India  (Southern  Asia) 

China  (Eastern  Asia) 
Australia  (Australia  &  New  Zealand) 
Peru  (South  America) 


MacDill  Air  Force  Base  (Florida,  USA) 
Osan  Air  Base  (South  Korea) 
McChord  Air  Force  Base  (Washington,  USA) 
Ramstein  Air-  Base  (Germany) 

Incirlik  Air  Base  (Turkey) 

Yokota  Air  Base  (Japan) 

Hickam  Air  Force  Base  (Hawaii,  USA) 
Charlotte  ANG  Base  (North  Carolina,  USA) 


Fig.  6  Optimal  military  HA/DR  aerial  delivery  supply  chain  network  design 


In  addition  to  this  solution  determined  by  incorporating  the  sb  probabilistic  scenario- 
iteration  combinations,  we  leveraged  the  Sample  Average  Approximation  (SAA) 
method29  within  the  MCDA  framework  to  better  estimate  the  optimal  solution  to 
the  stochastic  goal  program,  given  that  many  of  the  parameters  were  computed  by 
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randomly  sampling  from  a  prior  probability  distribution.  In  particular,  the  decision 
maker  used  the  MCDA  framework  to  run  the  optimization  model  for  30  separate 
instances  to  obtain  30  optimal  solutions;  in  theory,  the  average  of  these  solutions 
provides  an  unbiased  estimate  of  the  true  optimal  solution.  Given  that  each  of  the 
30  optimal  solutions  provided  different  results  (in  terms  of  which  storage  facilities 
to  open,  the  amount  of  inventory  to  store  at  each  depot  and  facility,  the  number  of 
trips  for  each  of  the  aerial  delivery  assets  pre-positioned,  and  the  amount  of  con¬ 
sumable  aid  to  procure  from  vendors),  it  did  not  make  practical  sense  to  provide 
the  decision  maker  with  mere  averages  of  the  decision  variable  values.  Instead,  the 
MCDA  framework  uses  the  results  of  the  SAA  method  to  construct  statistical  lower 
and  upper  bounds.30  This  approach  served  as  strategic  decision  support  by  help¬ 
ing  the  decision  maker  to  better  understand  the  range  of  possible  best-compromise 
solutions  to  the  multicriteria  military  HA/DR  aerial  delivery  supply  chain  network 
design  problem. 

For  example,  we  computed  statistical  lower  bounds  (LBs)  and  upper  bounds  (UBs) 
by  constructing  a  95%  confidence  interval  (Cl)  on  the  amount  of  consumable  aid 
inventory  to  store  at  each  of  the  depots  and  storage  facilities.  Table  7  shows  these 
results. 


Table  7  Statistical  upper  and  lower  bounds  for  amount  of  inventory  to  store 


Storage  Facilty/Depot  Name  95%  UB  Inventory 

95  %  LB  Inventory 

Ramstein  Air  Base  (Germany) 

634,227 

512,279 

Charlotte  ANG  Base  (North  Carolina,  USA) 

906,438 

801,936 

McChord  Air  Force  Base  (Washington,  USA) 

905,340 

777,940 

Hickam  Air  Force  Base  (Hawaii,  USA) 

770,074 

713,660 

Yokota  Air  Base  (Japan) 

588,970 

504,684 

Incirlik  Air  Base  (Turkey) 

642,841 

541,585 

Osan  Air  Base  (South  Korea) 

596,881 

506,079 

Ali  A1  Salem  Air  Base  (Kuwait) 

213,674 

71,340 

Soto  Cano  Air  Base  (Honduras) 

5.076 

-1,742 

MacDill  Air  Force  Base  (Florida,  USA) 

336,653 

177,053 

MCLB  Albany  (Georgia,  USA) 

MCLB  Base  Bars  tow  (California,  USA) 

7,419,858 

12,184,000 

3,750,276 

8,210,666 

The  decision  maker  used  the  results  in  Table  7  to  conduct  post-hoc  analyses  based 
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on  the  ranges  of  inventory  storage  amounts.  For  example,  the  95%  Cl  for  the  inven¬ 
tory  to  store  at  the  Soto  Cano  Air  Base  included  negative  values,  which  suggests 
that  this  storage  facility  should  probably  never  be  opened  for  pre -positioning  of  as¬ 
sets  and  inventory.  This  indicated  to  the  decision-maker  that  it  should  be  removed 
from  the  list  of  potential  candidate  storage  facilities. 

We  also  computed  the  statistical  LBs  and  UBs  on  the  amount  of  consumable  aid 
procured  from  the  vendors  for  storage  at  each  of  the  depots.  Table  8  provides  these 
results. 

The  decision  maker  used  the  results  in  Table  8  to  conduct  post-hoc  analyses  based 
on  the  ranges  of  the  inventory  procurement  amounts.  For  example,  the  inventory 
procurement  LB  for  vendors  VI  and  V2  were  significantly  lower  than  the  LB  for 
vendor  V3  for  inventory  storage  at  MCLB  Albany  (Georgia,  USA).  Also,  the  inven¬ 
tory  procurement  LB  for  vendor  V3  was  significantly  lower  than  the  LB  for  vendors 
VI  and  V2  for  inventory  storage  at  MCLB  Barstow  (California,  USA).  This  sug¬ 
gests  to  the  decision  maker  that  Albany  (Georgia,  USA)  may  want  to  consider  a 
single  sourcing  supplier  strategy  with  vendor  V3,  while  Barstow  (California,  USA) 
may  want  to  consider  a  multiple  sourcing  supplier  strategy  with  vendors  VI  and  V2. 


Table  8  Statistical  upper  and  lower  bounds  for  amount  of  inventory  to  procure 


Depot  Name 

Vendor 

95  %  UB  Procure 

95  %  LB  Procure 

MCLB  Albany  (Georgia,  USA) 

VI 

1,281,224 

314,510 

MCLB  Albany  (Georgia,  USA) 

V2 

2,376,233 

593,367 

MCLB  Albany  (Georgia,  USA) 

V3 

3,775,165 

2,829,635 

MCLB  Barstow  (California,  USA) 

VI 

5,020,958 

3,904,908 

MCLB  Barstow  (California,  USA) 

V2 

4,720,729 

2,831,271 

MCLB  Barstow  (California,  USA) 

V3 

2,449,673 

1,467,127 

We  also  computed  the  statistical  bounds  of  the  response  time,  cost,  and  demand 
unmet  to  further  explore  the  trade-offs  of  supply  chain  efficiency  versus  respon¬ 
siveness  (these  results  are  based  on  the  given  inputs  to  the  optimization  model).  In 
particular,  the  95%  Cl  of  the  average  response  time  (across  all  disaster  areas)  was 
[6,  7.21]  days,  the  average  total  supply  chain  cost  (across  all  stages  of  the  supply 
chain  network)  was  [$110,  $166]  billion,  and  the  average  total  demand  shortage 
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(across  all  disaster  areas)  was  [47.1,  47.2]  million  pounds.  These  results  provided  a 
baseline  for  "what  if"  sensitivity  analyses  conducted  for  the  decision  maker. 


4.4  Sensitivity  Analyses 

As  the  final  component  of  the  computational  experiment,  we  conducted  sensitivity 
analyses  to  investigate  how  changes  in  certain  model  parameters  affected  the  opti¬ 
mal  solution  as  a  means  to  provide  strategic  decision  support  of  the  supply  chain 
network  trade-offs. 

We  first  adjusted  the  data  parameter  R,  which  is  the  decision-maker’s  target  total 
flight  time  budget  that  directly  affects  the  maximum  number  of  trips  with  consum¬ 
able  aid  from  storage  facility  j  to  disaster  area  i  via  aerial  delivery  mode  k  (the 
maximum  number  of  trips  increases  as  R  increases).  Here,  all  other  model  parame¬ 
ters  were  fixed  as  previously  described,  but  we  evaluated  the  changes  in  the  average 
response  time,  supply  chain  cost,  and  demand  shortage  as  R  was  varied  from  1  day 
(24  h)  to  14  days  (336  h).  We  first  noted  that  the  optimal  subset  of  storage  facilities 
to  open  changed  at  R  =  4  days  (96  h),  where  the  optimization  model  decided  to 
open  the  Soto  Cano  Air  Base  (Honduras)  for  pre-positioning  while  closing  the  Osan 
Air  Base  (South  Korea). 

In  addition,  we  found  that  as  the  value  R  increased,  the  average  response  time 
(across  all  disaster  areas)  tended  to  decrease.  The  average  total  supply  chain  cost, 
however,  tended  to  increase  through  R  =  7  days  but  then  decreased  slightly.  More¬ 
over,  the  average  total  demand  shortage  tended  to  decrease  through  R  —  7  days  but 
then  increased  slightly — a  relationship  inverse  to  that  of  total  cost.  Table  9  displays 
these  results. 
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Table  9  Supply  chain  network  design  trade-offs  from  adjusting  R 


R 

(d) 

R 

(h) 

Time 

(d) 

Total  Cost 

(000s) 

Demand  Unmet 

(lb) 

1 

24 

6.04 

$  65,699,925 

47,518,280 

2 

48 

6.08 

$  134,834,650 

47,223,830 

3 

72 

6.11 

$  152,475,550 

47,124,980 

4 

96 

6.01 

$  157,934,050 

47,093,180 

5 

120 

5.48 

$  158,980,950 

47,083,630 

6 

144 

5.35 

$  172,945,300 

47,031,580 

7 

168 

5.21 

$  174,846,850 

47,012,530 

8 

192 

5.18 

$  173,146,700 

47,018,930 

9 

216 

4.92 

$  150,118,800 

47,134,130 

10 

240 

4.92 

$  143,208,850 

47,166,130 

11 

264 

4.93 

$  141,531,750 

47,190,130 

12 

288 

4.82 

$  130,832,750 

47,233,180 

13 

312 

4.92 

$  178,315,600 

47,000,730 

14 

336 

4.85 

$  136,600,800 

47,202,780 

The  scatter  plots  in  Fig.  7  also  depicts  these  supply  chain  network  design  trade-offs 
of  average  response  time,  cost,  and  demand  unmet.  In  each  of  the  plots,  we  also 
show  the  LOWESS  smoother  that  uses  locally  weighted  polynomial  regression  to 
see  the  overall  trend. 

The  scatter  plots  indicate  the  direct  trade-off  between  cost  and  demand  unmet.  For 
example,  when  R  =  7  days  the  average  total  cost  reached  its  maximum  while 
the  average  amount  of  demand  unmet  reached  its  minimum.  These  plots  reveal  an 
inverse  relationship  between  supply  chain  efficiency  and  responsiveness.  Note  that 
at  R  =  7  days,  the  average  total  response  time  is  not  at  its  minimum  or  maximum 
point. 
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Fig.  7  Scatter  plots  of  supply  chain  network  design  trade-offs  from  adjusting  R 


To  further  illustrate  the  supply  chain  network  design  trade-offs,  Fig.  8  depicts  the 
bicriteria  objective  spaces  between  time  and  cost,  time  and  demand,  and  cost  and 
demand;  these  plots  highlight  the  nondominated  solution  sets. 

For  example,  in  the  left  plot  of  Fig.  8  there  are  2  nondominated  solutions  (4.82  d, 
$130,832,750,000)  at  R  =  12  days  and  (6.04  d,  $65,699,925,000)  at  R  =  1  day. 
The  nondominated  solutions  are  also  evident  in  the  middle  plot.  In  the  right  plot, 
there  appears  to  be  a  linear  trade-off  between  average  demand  unmet  and  average 
total  supply  chain  cost;  this  confirms  our  previous  finding  that  these  supply  chain 
network  design  objectives  are  inversely  related. 

Next,  we  adjusted  the  data  parameter  T,  which  is  the  decision-maker’s  target  total 
flight  time  budget  that  directly  affects  the  maximum  number  of  trips  with  consum¬ 
able  aid  from  depot  d  to  storage  facility  j  via  the  C-17  aerial  delivery  platform; 
again,  note  that  the  maximum  number  of  trips  increases  as  T  increases.  Here,  all 
other  model  parameters  were  fixed  as  previously  described  in  the  setup  of  the  com¬ 
putational  experiment,  but  we  evaluated  the  changes  in  the  average  response  time, 
supply  chain  cost,  and  demand  shortage  as  T  was  varied  from  2  days  (48  h)  to  14 
days  (336  h).  We  first  noted  that  the  optimal  subset  of  storage  facilities  to  open 
changed  several  times  at  T  =  2,3,4  days;  the  storage  facilities  that  changed  in¬ 
cluded  Soto  Cano  (Honduras),  Osan  (South  Korea),  and  Ali  A1  Salem  (Kuwait)  Air 
Bases. 
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Average  Respone  Time  (Days) 


Average  Response  Time  (Days) 


Average  Demand  Unmet  (lbs) 


Fig.  8  Supply  chain  network  design  bi-criteria  solution  sets  from  adjusting  R 


Moreover,  we  found  that  as  the  value  T  increased,  both  the  average  response  time 
and  the  average  total  supply  chain  cost  tended  to  increase;  however,  the  amount  of 
demand  unmet  tended  to  decrease.  Again,  there  appears  to  be  an  inverse  relationship 
between  cost  and  demand  unmet.  Table  10  displays  these  results. 
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Table  10  Supply  chain  network  design  trade-offs  from  adjusting  T 


T 

(d) 

T 

(h) 

Time 

(d) 

Total  Cost 

(000s) 

Demand  Unmet 

(lb) 

2 

48 

5.33 

$56,041,825 

47,579,104 

3 

72 

5.66 

$91,068,740 

47,409,980 

4 

96 

5.73 

$  110,316,950 

47,321,980 

5 

120 

6.11 

$  152,475,550 

47,124,980 

6 

144 

6.05 

$  184,322,650 

46,973,430 

7 

168 

5.87 

$  194,610,500 

46,927,180 

8 

192 

5.87 

$201,143,350 

46,896,180 

9 

216 

5.86 

$  195,205,400 

46,909,530 

10 

240 

5.94 

$  198,269,300 

46,894,230 

11 

264 

5.96 

$  208,352,800 

46,865,030 

12 

288 

5.97 

$  209,670,200 

46,861,280 

13 

312 

5.99 

$  209,277,500 

46,861,280 

14 

336 

5.99 

$209,041,400 

46,861,280 

The  scatter  plots  with  LOWESS  smoothers  in  Fig.  9  further  illustrate  these  supply 
chain  network  design  trade-offs. 

In  the  left  plot,  the  average  response  time  peaks  at  T  =  5  days,  indicating  to  the 
decision  maker  that  the  target  total  flight  time  budget  should  be  somewhere  between 
2  to  4  days  or  7  to  14  days  (depending  on  the  desired  trade-offs).  In  the  middle  and 
right  plots,  the  distinct  trade-off  between  cost  and  demand  unmet  is  again  clearly 
illustrated.  Figure  10  depicts  the  objective  spaces  between  these  criteria. 
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Fig.  9  Scatter  plots  of  supply  chain  network  design  trade-offs  from  adjusting  T 


Average  Respone  Time  (Days) 


Average  Response  Time  (Days) 


Average  Demand  Unmet  (lbs) 


Fig.  10  Supply  chain  network  design  bi-criteria  solution  sets  from  adjusting  T 


The  left  plot  in  Fig.  10  illustrates  that  average  response  time  and  total  cost  are  not 
conflicting,  as  there  is  only  one  nondominated  solution  (5.33  d,  $56,041,825,000) 
at  T  =  2  days.  On  the  other  hand,  the  middle  plot  contains  numerous  nondominated 
and  dominated  solutions  between  the  time  and  demand  unmet  criteria;  for  example, 
the  bicriteria  solution  (5.33  d,  47,579,104  lb)  at  T  =  2  days  is  nondominated  while 
(6.1 1  d,  47,124,980  lb)  at  T  =  5  days  is  dominated.  The  right  plot  again  shows  the 
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almost  perfectly  linear  trade-off  between  average  demand  unmet  and  average  total 
supply  chain  cost. 

Given  that  the  MCDA  framework  employed  weighted  goal  programming,  there  was 
tremendous  sensitivity  in  the  assignment  of  decision-maker  weights  for  each  goal. 
Because  the  first  2  goals  wanted  to  minimize  response  time  and  minimize  cost, 
an  optimal  solution  would  not  open  any  storage  facilities.  On  the  other  hand,  the 
third  goal  wanted  to  maximize  the  amount  of  demand  met  (or  minimize  amount  of 
demand  unmet);  thus,  an  optimal  solution  would  open  as  many  storage  facilities  as 
possible.  Table  1 1  shows  the  supply  chain  network  design  trade-offs  (response  time, 
supply  chain  cost,  and  demand  unmet)  given  different  combinations  of  decision¬ 
maker  weights. 


Table  11  Supply  chain  network  design  trade-offs  from  varying  decision-maker  weights 


W1 

W2 

W3 

W4 

Time 

(d) 

Total  Cost 

(000s) 

Demand  Unmet 

(lb) 

1 

1 

too 

iOO.OOO 

6.11 

$  152,475,550 

47,124,980 

1 

2 

too 

fOO.OOO 

6.30 

$  89,833,355 

47,406,130 

f,000 

1 

1 

iOO.OOO 

6.06 

$  145,540,950 

47,161,780 

too 

1 

1 00,000 

i.OOO 

6.04 

$  148,996,700 

47,133,330 

From  Table  11,  we  found  that  the  third  goal  required  significantly  greater  weights 
than  the  first  and  second  goals  to  achieve  a  pragmatic  best  compromise  solution. 
However,  the  weight  of  the  second  goal  was  the  most  significant  in  its  effect  on  the 
solution.  For  example,  when  increasing  the  weight  of  the  second  goal  by  merely 
one,  the  average  response  time  increased  by  0.19  days,  total  cost  decreased  by  over 
$62  billion,  and  amount  of  demand  unmet  increased  by  over  28,000  lb.  As  another 
instance,  when  increasing  the  weight  of  the  first  goal  by  999  (and  decreasing  the 
weight  of  the  third  goal  by  99),  the  average  total  aerial  delivery  response  time  de¬ 
creased  by  0.05  days,  total  cost  decreased  by  nearly  $7  billion,  and  demand  unmet 
increased  by  over  36,000  lb.  As  a  result  of  these  findings,  the  selection  of  appropri¬ 
ate  goal  weights  by  the  decision  maker  requires  significant  attention. 
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5.  Concluding  Remarks 


In  this  report,  we  developed  a  multiple  criteria  decision  analysis  framework  to  sup¬ 
port  the  strategic  logistics  planning  and  cost  assessment  of  military  HA/DR  aerial 
delivery  operations.  The  decision  support  tool  facilitated  both  responsive  and  effi¬ 
cient  logistics  planning  under  uncertainty  by  optimizing  the  military  HA/DR  aerial 
delivery  supply  chain  network. 

In  particular,  the  framework  employed  stochastic,  mixed-integer,  weighted  goal 
programming  to  optimize  supply  chain  network  design,  logistics  costs,  staging  lo¬ 
cations,  procurement  amounts,  and  inventory  levels.  Moreover,  the  framework  in¬ 
corporated  several  methodologies  to  provide  robustness  and  mitigate  uncertainty  in 
the  system.  Finally,  the  framework  enabled  the  decision  maker  to  explore  the  supply 
chain  network  design  trade-offs  of  average  total  aerial  delivery  response  time,  total 
supply  chain  cost,  and  amount  of  demand  unmet. 

We  determined  a  recommended  network  design  based  on  the  criteria  and  data  pro¬ 
vided  by  NSRDEC,  which  provided  an  average  response  within  6  days  and  a  supply 
chain  cost  of  $153  billion.  Sensitivity  analysis  on  the  solution  showed  a  direct  trade¬ 
off  between  cost  and  unmet  demand.  Further,  as  the  flight  time  budget  increased, 
the  average  response  time  and  supply  chain  cost  increased  while  unmet  demand 
decreased.  The  use  of  goal  programming,  however,  gave  significant  flexibility  for 
the  decision  maker.  Although  we  applied  the  modeling  framework  to  a  particular 
application,  it  can  easily  be  adopted  to  other  HA/DR  aerial  delivery  applications. 

There  are  several  limitations  to  the  work  that  should  be  mentioned.  First,  parameter 
estimation  for  rare  events  is  difficult  since  the  historical  data  tends  to  be  sparse.  In 
addition,  we  lacked  information  to  be  able  to  define  detailed  distributions  on  several 
factors,  including  vendor  procurement  costs  and  cruise  speed,  and  instead  modeled 
a  range  of  values  using  the  uniform  distribution.  We  used  average  values  derived 
by  the  regression  model  as  the  basis  for  demand  the  scenarios.  The  results  might 
change  if  we  had  data  to  better  model  potential  extremes.  In  addition,  the  regres¬ 
sion  model  used  to  estimate  the  affected  population  per  year  only  explained  21%  of 
the  variation,  in  part  due  to  the  limited  number  of  factors  in  the  EM-DAT  data  set. 
Further,  the  commodities  modeled  here  are  simple  in  that  they  are  not  perishable 
and  may  be  characterized  purely  by  their  weight.  Other  applications  may  require 
more  sophisticated  models  of  capacity  that  take  into  account  volume,  weight,  and 
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perishability.  Finally,  we  minimized  the  sum  of  the  expected  goal  deviations  in  the 
optimization  model.  Although  this  is  a  very  common  form  in  practice  and  in  previ¬ 
ous  literature,  using  expectations  has  limitations  when  the  variance  is  quite  high.31 
It  also  assumes  that  the  decision  makers  are  risk  neutral. 

In  future  work,  we  plan  to  conduct  additional  sensitivity  analyses  by  adjusting  the 
levels  of  the  experimental  design  factors.  We  also  intend  to  evaluate  the  effects  of 
procurement  cost  variability,  as  the  range  set  by  the  decision  maker  for  this  study 
may  have  been  artificially  tight.  Further,  we  aim  to  reconfigure  the  optimization 
model  formulation  by  incorporating  a  fourth  decision-maker  goal  that  accounts  for 
the  conflicting  criterion  of  supply  chain  disruption/risk  as  well  as  altering  some 
of  the  model  constraints  (such  as  allowing  storage  facilities  to  serve  more  than 
one  disaster  area).  Finally,  we  plan  on  tying  the  strategic  model  developed  here  to 
operational  planning  models. 
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Appendix  A.  Screenshots  of  MCDA  Framework 
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Military  HA/DR  Aerial  Delivery  Cost  Assessment  and  Logistics  Modeling  Tool 

Multiple  Criteria  Decision  Analysis  (MCDA)  Framework 


Write  Input  Files 


File  Description 

This  file  provides  a  data  interface  between  Excel  and  the  General  Algebraic  Modeling  System  (GAMS)  using  VBA. 

Write  Button 

-  Writes  out  input  data  to  read  into  GAMS. 


Optimize 


Optimize  for  Statistical  Bounds 

Optimize  Button 

-  Opens  a  Command  Window  and  invokes  the  'Run.bat'  batch  file  to  execute  GAMS.  Optimize  over  30  Model 

-  GAMS  then  solves  the  optimization  model  (one  instance)  and  outputs  csv  data  files.  Instances 


Read  Results  Files 


Read  Button 

-  Reads  in  GAMS  model  output  and  conducts  statistical  analysis. 


pennState 


The  Marcus  Department  of 
Industrial  and  Manufacturing 
Engineering 


MliL 


Dash  Board 


22^21  Model  Output  Charts  Statistical  Bounds 


Sensitivity  Analysis 
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Fig.  A-l  MCDA  tool  dashboard 
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Candidate  Storage  Facilities 

Degrees 

Radians 

Storage  Facility 

Name 

Location 

Latitude 

Longitude 

Latitude 

Longitude 

"FI" 

Ramstein  Air  Base 

Germany 

49443889 

7.602222 

0  86295866  0.132683804 

"F2" 

Charlotte  Air  National  Guard  Base 

North  Carolina,  USA 

35.216111 

-80.931944 

0.61463709 

-1  41252889 

"F3" 

McChord  Air  Force  Base 

Washington.  USA 

47  137778 

-122.47639 

0  82270943  -2  13761624 

"F4" 

Hickam  Air  Force  Base 

Hawaii.  USA 

21  318611 

-157.9225 

0  37207995 

-2  75626759 

"F5" 

Yokota  Air  Base 

Japan 

35  748611 

139.348611 

0  62393096  2  43209207 

"F6" 

Incirlik  Air  Base 

Turkey 

37.001944 

35.425833 

0  64580575  0.618297426 

"FT 

Osan  Air  Base 

South  Korea 

37  090556 

127  029722 

0  64735232  2.217086897 

”F8" 

Ali  Al  Salem  Air  Base 

Kuwait 

29  346667 

47.520556 

0  51219596  0  829390165 

"F9" 

Soto  Cano  Air  Base 

Honduras 

14.3825 

-87.621111 

0.25102198 

-1  52927688 

"F10" 

MacDill  Air  Force  Base 

Florida.  USA 

27.849444 

-82.521111 

0.48606449 

-1.44026509 

Depots 

Degrees 

Radians 

Depot 

Name 

Shipping  Location 

Latitude  Longitude 

Latitude  Longitude 

“D1" 

”D2" 

Marine  Corps  Logistics  Base  Albany 
Marine  Corps  Logistics  Base  Barstow 

Southwest  Georgia  Regional  Airport  (ABY),  Georgia.  USA 
Los  Angeles  International  Airport  (LAX).  California,  USA 

31.535556  -84.194444 
33.9425  -118.40806 

0.55039928  -1  46947026 
0.59240838  -2  06661044 

Potential  Disaster  Areas 

Deg 

rees 

Radians 

Disaster  Area 

Name  (Geographic  Region) 

Location  of  Largest  City 

Latitude 

Longitude 

Latitude 

Longitude 

"Al" 

Haiti  (Caribbean) 

Port-au-Prince 

18.533333 

-72.333333 

0 32346768 

-1  26245482 

"A2" 

Indonesia  (South-Eastern  Asia) 

Jakarta 

-6.175 

106.828333 

-0.1077741 

1  864506145 

"A3" 

Mexico  (Central  America) 

Mexico  City 

19.433333 

-99.133333 

0  33917565 

-1.73020306 

"A4" 

Tanzania  (Eastern  Africa) 

Dar  es  Salaam 

-6.8 

39.283333 

-0.1186824 

0.685623502 

"A5" 

India  (Southern  Asia) 

Mumbai 

18.975 

72.825833 

0.33117623 

1.271050566 

"A6" 

China  (Eastern  Asia) 

Shanghai 

31.2 

121.5 

0.54454273 

2.120575041 

"A7" 

Australia  (Australia  &  New  Zealand) 

Sydney 

-33.85997 

151  209444 

-0.590968 

2  639102658 

"A8" 

Peru  (South  America) 

Lima 

-12  04333 

-77  028333 

-0  2101958 

-1  34439803 

Instructions: 

1)  You  may  adjust  the  name,  location,  and  latitude/longitude  (degrees) 
of  the  candidate  storage  facilities. 

2)  You  may  adjust  the  name,  shipping  location,  and  latitude/longitude 
(degrees)  of  the  depots.. 

3)  You  may  change  the  number  of  aerial  delivery  modes  positioned  at 
each  candidate  storage  facility. 


#  of  Aerial  Delivery  Modes  Positioned  at  Facility 

FI 

C17 

2 

C130 

3 

CH47 

5 

F2 

3 

2 

0 

F3 

2 

3 

3 

F4 

3 

3 

4 

F5 

2 

2 

2 

F6 

2 

2 

0 

F7 

2 

3 

3 

F8 

1 

1 

0 

F9 

0 

0 

5 

F10 

1 

1 

0 

Dash  Board  Data  Input  Model  Output  Charts  Statistical  Bounds 


Sensitivity  Analysis 
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Fig.  A-2  MCDA  tool  data  input 


Opened  (1=Yes,  0=No) 

Facility  Name 

Inventory 

Number  of  Trips  by  Aerial  Delivery  Mode 

1 

Ramstein  Air  Base 

665600 

4 

2 

0 

1 

Charlotte  Air  National  Guard  Base 

832000 

3 

5 

0 

1 

McChord  Air  Force  Base 

896000 

4 

5 

0 

1 

Hickam  Air  Force  Base 

769600 

3 

2 

1 

1 

Yokota  Air  Base 

640000 

3 

5 

0 

1 

Incirlik  Air  Base 

588800 

3 

4 

0 

1 

Osan  Air  Base 

563200 

2 

4 

0 

0 

Ali  Al  Salem  Air  Base 

0 

0 

0 

0 

0 

Soto  Cano  Air  Base 

0 

0 

0 

0 

1 

MacDill  Air  Force  Base 

435200 

4 

7 

0 

"Cl  7" 

"C130" 

"CH47" 

Depot  Name 

Inventory 

Marine  Corps  Logistics  Base  Albany 

14208000 

Depot 

Facility 

Number  of  Trips  by  C17 

Marine  Corps  Loqistics  Base  Barstow 

2496000 

”D1" 

"FI" 

6 

"Dr 

"F2" 

8 

Disaster  Area  (Geographic  Region)  Name 

Servicing  Storage  Facility 

"Dr 

"F3" 

7 

Haiti  (Caribbean) 

"F10” 

"D1" 

"F4" 

6 

Indonesia  (South-Eastern  Asia) 

"F7" 

„Di” 

"F5" 

5 

Mexico  (Central  America) 

"F3" 

-or 

"F6" 

5 

Tanzania  (Eastern  Africa) 

"FI" 

"D1” 

■FT” 

4 

India  (Southern  Asia) 

"F6" 

"D1” 

"F8" 

0 

China  (Eastern  Asia) 

"F5" 

"D 1 " 

"F9" 

0 

Australia  (Australia  &  New  Zealand) 

"F4" 

"D1" 

"F10" 

7 

Peru  (South  America) 

”F2" 

"D2" 

"FI” 

5 

"D2" 

"F2" 

5 

Vendor 

Amount  Procured 

"D2" 

"F3” 

7 

■vr 

3072000 

"D2" 

"F4" 

7 

”V2" 

5568000 

„D2" 

"F5" 

5 

”V3" 

5568000 

"D2" 

"F6" 

5 

"VI” 

2496000 

"D2" 

„F7" 

5 

"V2" 

0 

"D2" 

"F8" 

0 

“V3” 

0 

”D2" 

"F9" 

0 

"D2" 

"F10" 

0 

Goal 

Average  Under 

Average  Over 

■T 

0.00 

21715.88 

"2” 

0.00 

152475550000.00 

”3" 

47124979.95 

4400662.43 

"4" 

0.00 

47124979.95 

Average  Response  Time  Across  All  Disaster  Areas 

Hours  =  146  73 
Days  =  6.11 

Average  Total  Supply  Chain  Cost  in  the  Full  Network 

$  152,475,550,000.00 


Average  Total  Demand  Shortage  Across  All  Disaster  Areas 

Pounds  =47124979  95 


Goal 


"2" 

•T 


Median  Under 

0.00 
0.00 
0.00 
0  00 


Median  Over 


21715.88 
152475550000.00 
5234551  70 
000 


Median  Response  Time  Across  All  Disaster  Areas 

Hours  =  146.73 
Days  =  6.11 

Median  Total  Supply  Chain  Cost  in  the  Full  Network 

$  152,475,550,000.00 

Median  Total  Demand  Shortage  Across  All  Disaster  Areas 

Pounds  =  0.00 

Percent 

#  of  Times  Demand  Unmet:  8  13% 
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Fig.  A-3  MCDA  tool  model  output 
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Number  of  Trips  by  Aerial  Delivery  Mode  to  Disaster  Areas 

8 


Amount  of  Inventory  to  Store  at  Depot  or  Facility 


T>2" 

■or 


"F8" 

“FT' 

”F6” 

”F5" 

"F4" 

"F3" 

•FT 

"FI" 


4000000 


6000000 


10000000  12000000  14000000  16000000 


"FI"  TT  ”F3"  T4"  T5"  "F6"  T  T  "F8"  "F9"  "F10"  "D1"  "D2" 

■  Inventory  665600  832000  896000  769600  640000  588800  563200  0  0  435200  14208000  2496000 


Number  of  Trips  by  C17  from  "D1"  to 
Storage  Facilitites 


"FI"  "F2"  T3"  "F4"  "F5"  "F6"  T7"  "F8"  T9"  "F10" 


■  Trips  6876554007 

Number  of  Trips  by  C17  from  "D2”  to 
Storage  Facilitites 


Dash  Board  Model  Output  ]  Charts  Statistical  Bounds 


Sensitivity  Analysis 
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Fig.  A-4  MCDA  tool  output  charts 


146706400000.00 

4424065.76 

47121636.62 


4713644529 


95%  LB  Pr 


138-0C2  834998 
3655340  1541460 
4557965  3429235 
4722206  3572194 
4010882  1902718 
2118943 _ 1004257 


Fig.  A-5  MCDA  tool  statistical  bounds 
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R  (Days)  R  (Hours)  Time  (Days)  Total  Cost  (000s)  Demand  Unmet  (lbs) Remarks 


1 

24 

6.04 

$ 

65,699,925 

47518280 

2 

48 

6.08 

S 

134.834,650 

47223830 

3 

72 

6.11 

$ 

152,475,550 

47124980 

4 

96 

6.01 

s 

157,934,050 

470931 80  Soto  Cano  Opened  but  Osan  Closed 

5 

120 

548 

s 

158.980,950 

47083630 

6 

144 

5.35 

$ 

172,945,300 

47031580 

7 

168 

5.21 

s 

174.846,850 

47012530 

8 

192 

5.18 

$ 

173,146,700 

47018930 

9 

216 

4.92 

$ 

150,118,800 

47134130 

10 

240 

4.92 

s 

143,208,850 

47166130 

11 

264 

493 

s 

141,531,750 

47190130 

12 

288 

4.82 

s 

130,832,750 

47233180 

13 

312 

4.92 

s 

178,315,600 

47000730 

14 

336 

4.85 

$ 

136,600,800 

47202780 

T  (Days)  T  (Hours)  Time  (Days) 

Total  Cost  (000s) 

Demand  Unmet  (lbs)  Remarks 

2 

48 

5.33 

$ 

56,041,825 

47579104  Soto  Cano  Opened  but  Osan  Closed 

3 

72 

5.66 

$ 

91,068,740 

47409980  Osan  Opened  but  Ali  Al  Salem  Closed 

4 

96 

5.73 

S 

110,316,950 

47321980  Ali  Al  Salem  Opened  by  Osan  Closed 

5 

120 

6.11 

s 

152,475,550 

47124980 

6 

144 

6.05 

$ 

184,322,650 

46973430 

7 

168 

5.87 

s 

194.610,500 

46927180 

8 

192 

5.87 

$ 

201,143,350 

46896180 

9 

216 

5.86 

$ 

195,205,400 

46909530 

10 

240 

5.94 

$ 

198,269.300 

46894230 

11 

264 

596 

s 

208,352,800 

46865030 

12 

288 

5.97 

s 

209.670,200 

46861280 

13 

312 

5.99 

s 

209,277,500 

46861280 

14 

336 

5.99 

s 

209,041,400 

46861280 

W1  W2 

W3 

W4 

Time  (Days)  Total  Cost  (000s) 

Demand  Unmet  (lbs) 

1  1 

100 

100000 

6.11 

S 

152,475,550 

47124980 

1  2 

100 

100000 

6.30 

s 

89,833,355 

47406130 

1000  1 

1 

100000 

6.06 

s 

145,540,950 

47161780 

100  1 

100000 

1000 

6.04 

s 

148,996,700 

47133330 

Fig.  A-6  MCDA  tool  sensitivity  analysis 
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Intentionally  left  blank. 
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Appendix  B.  GAMS  Programming  Code  -  Single  Model  Instance 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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$Title  HA/DR  Aerial  Delivery  Optimization  Model 
$ontext 

This  multiple  criteria  decision  analysis  (MOD A)  framework  is  used 
for  humanitarian  assistance  /  disaster  relief  aerial  delivery  logistics 
planning  and  cost  assessment. 

$offtext 

*Set  pseudo-random  number  generator  seed 
execseed  =  le8*(frac(jnow)); 

*  Limits  the  amount  of  CPLEX  solver  time  to  100,000  seconds 
Option  ResLim=  100000; 

*  Limits  the  number  of  CPLEX  solver  iterations  to  100,000  runs 
Option  IterLim=  100000; 

*  Limits  the  number  of  decimal  places  to  4 
Option  decimals=4; 


Sets 

v  consumable  aid  vendors  (suppliers)  /VI,  V2,  V3/ 
d  consumable  aid  depots  /D 1 ,  D2 / 

j  candidate  storage  facilities  for  pre-positioning  /FI,  F2,  F3,  F4,  F5,  F6,  F7, 
F8,  F9,  F 1 0/ 

i  potential  affected  disaster  areas  (demand  regions)  /Al,  A2,  A3,  A4,  A5,  A6, 
A7,  A8/ 

k  aerial  delivery  modes  /Cl 7,  Cl 30,  CH47/ 

g  goals /I  *4/ 

s  demand  scenarios  / 1  *  8/ 

me  trials  for  Monte  Carlo  Simulation  /I  *  1 000/ 

fact  DOE  factors  generating  each  scenario  /prob,  dislmp,  DaysFood, 

SCDisrupt/ 

b  model  iterations  / 1  *  8/; 

Parameters 

w(g)  decision-maker  weight  for  goal  g 

Dem(i,s,b)  expected  consumable  aid  demand  (lbs)  per  day  at  disaster  area  i  in 
scenario  s 

ds(i,j)  geodesic  distance  (NM)  from  storage  facility  j  to  disaster  area  i 

l(d,j)  geodesic  distance  (NM)  from  depot  d  to  storage  facility  j 

Spd(k)  average  speed  (knots)  of  aerial  delivery  mode  k 
LoadDelay(k)  load  time  for  aid  (per  trip)  of  aerial  delivery  mode  k 
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t(i,j,k)  travel  time  (per  trip)  from  facility  j  to  disaster  area  i  via  delivery  mode 
k 

r(d,j)  travel  time  (per  trip)  from  depot  d  to  storage  facility  j 

OMF(k)  operation  &  maintenance  &  fuel  cost  ($  per  hr)  of  aerial  delivery 

mode  k 

TDcost(d,j)  transportation  cost  (per  trip)  from  depot  d  to  storage  facility  j 
TScost(i,j,k)  transportation  cost  (per  trip)  from  facility  j  to  disaster  area  i  via 
delivery  mode  k 

Atd(d,j,s,b)  amount  of  aid  (lbs)  transported  (per  trip)  from  depot  d  to  storage 
facility  j  in  scenario  s  and  iteration  b 

Ats(i,j,k,s,b)  amount  of  aid  (lbs)  transported  (per  trip)  from  storage  facility  j  to 
disaster  area  i  via  delivery  mode  k  in  scenario  s  and  iteration  b 
IHDcost(d)  inventory  holding  cost  (per  lb)  at  depot  d 
IHScost(j)  inventory  holding  cost  (per  lb)  at  storage  facility  j 
Pcost(v)  procurement  cost  (per  lb)  of  consumable  aid  from  vendor  v 
Cap(j)  inventory  holding  capacity  (lbs)  at  storage  facility  j 

Cp(d)  inventory  holding  capacity  (lbs)  at  depot  d 

Fcost(j)  fixed  cost  of  opening  storage  facility  j 

Ntd(d,j)  maximum  number  of  trips  with  consumable  aid  from  depot  d  to 
storage  facility  j 

Nts(i,j,k)  maximum  number  of  trips  with  consumable  aid  from  storage  facility 
j  to  disaster  area  i  via  aerial  delivery  mode  k 

ADM(j,k)  number  of  aerial  delivery  mode  k  pre-positioned  at  storage  facility 

j 

its(b,  s,  fact)  DOE  factors  for  each  scenario  and  iteration  combination 
Di(i,s)  ddd 


*  Input  Data 


*  Decision-maker  weight  for  goal  g 
w(T)  =  1; 

w('2')  =  1; 
w(’3’)=  100; 
w('4')=  100000; 

*  Expected  consumable  aid  demand  (lbs)  per  day  at  disaster  area  i  in  scenario  s 
Table  Di(i,s) 

$ondelim 

$include  DemandDisaster.csv 
$offdelim 


*  Expected  consumable  aid  demand  (lbs)  per  day  at  disaster  area  i  in  scenario  s 
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loop(i, 

loop(s, 

loop(b, 

Dem(i,s,T)  =  Di(i,s); 

Dem(i,s,'2')  =  Di(i,s); 

Dem(i,s,'3')  =  Di(i,s); 

Dem(i,s,'4')  =  Di(i,s); 

Dem(i,s,'5')  =  Di(i,s); 

Dem(i,s,'6')  =  Di(i,s); 

Dem(i,s,'7')  =  Di(i,s); 

Dem(i,s,'8')  =  Di(i,s); 

););); 

*  Geodesic  distance  (NM)  from  storage  facility  j  to  disaster  area  i 
Table  ds(i,j) 

Sondelim 

$include  Dist_StageDisaster.csv 
Soffdelim 


*  Geodesic  distance  (NM)  from  depot  d  to  storage  facility  j 
Table  l(d,j) 

$ondelim 

$include  Dist_DepotStage.csv 
Soffdclim 


^Conduct  a  Monte  Carlo  Simulation  to  determine  Spd(k)  and  LoadDelay(k) 
Parameters 

count  1  counter  for  the  Monte  Carlo  Simulation 

count2  counter  for  the  Monte  Carlo  Simulation 

count3  counter  for  the  Monte  Carlo  Simulation 

count4  counter  for  the  Monte  Carlo  Simulation 

count5  counter  for  the  Monte  Carlo  Simulation 

count6  counter  for  the  Monte  Carlo  Simulation 

vel(k,mc)  velocity  computation  for  Monte  Carlo  Simulation 

delay(k,mc)  load  delay  time  computation  for  Monte  Carlo  Simulation 

count  1  =  0; 
count2  =  0; 
count3  =  0; 
count4  =  0; 
count5  =  0; 
count6  =  0; 
loop(k, 
loop(mc, 
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*  Uniform  distribution  from  400  to  500  knots  (cruise  =  450  knots)  -  C17 
vel('C17',mc)  =  uniform(400,500); 

count  1  =  count  1  +  vel('C17’,mc); 

*  Uniform  distribution  from  250  to  320  knots  (cruise  =  292  knots)  -  Cl 30 
vel('C130’,mc)  =  uniform(250,320); 

count2  =  count2  +  vel('C130’,mc); 

*  Uniform  distribution  from  90  to  170  knots  (cruise  =130  knots)  -  CH47 
vel('CH47',mc)  =  uniform(90,170); 

count3  =  count3  +  vel('CH47',mc); 

*  Uniform  distribution  from  1.75  to  3.25  hours  -  C17 
delay('C17’,mc)  =  uniform(1.75,3.25); 

count4  =  count4  +  delay(’C17’,mc); 

*  Uniform  distribution  from  0.75  to  2.25  hours  -  C130 
delay('C130’,mc)  =  unifonn(0.75,2.25); 

count5  =  count5  +  delay(’C130’,mc); 

*  Uniform  distribution  from  0.25  to  1.50  hours  -  CH47 
delay('CH47’,mc)  =  unifonn(0.25,1.5); 

count6  =  count6  +  delay('CH47’,mc); 

); 

*  Average  cruise  speed  (knots)  of  aerial  delivery  mode  k 
Spd('C17')  =  count  l/card(mc); 

Spd('C130’)  =  count2/card(mc); 

Spd('CH47')  =  count3/card(mc); 

*  Load  delay  time  (hr)  of  aerial  delivery  mode  k 
LoadDelay('C17')  =  count4/card(mc); 

LoadDelay('C130')  =  count5/card(mc); 

LoadDelay('CH47')  =  count6/card(mc); 

); 


*  Travel  time  (per  one-way  trip)  from  facility  j  to  disaster  area  i  via  delivery  mode 
k 

loop(i, 

loop(j, 

loop(k, 

t(i,j,k)  =  (ds(i,j)  /  Spd(k))  +  LoadDelay(k); 

); ); ); 

*  Travel  time  (per  one-way  trip)  from  depot  d  to  storage  facility  j 
loop(d, 

loop(j, 

r(d,j)  =  (l(d,j)  /  Spd(’C17'))  +  LoadDelay(’C17'); 

); ); 
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*  Operation,  Maintenance,  and  Fuel  Cost  in  $  per  hour 
OMF('C17')  =uniform(  12000, 13500); 

OMF('C130')  =  unifonn(4000,4700); 

OMF('CH47')  =  unifonn(9000, 10300); 

*  Transportation  cost  (per  trip)  from  depot  d  to  storage  facility  j 
loop(d, 

loopG, 

TDcost(d,j)  =  (250  +  OMF(’C17’))  *  r(d,j); 

); ); 

*  Transportation  cost  (per  trip)  from  facility  j  to  disaster  area  i  via  delivery  mode 
k 

loop(i, 

loopG, 

loop(k, 

TScost(i,j,k)  =  (250  +  OMF(k))  *  t(i,j,k); 

););); 


*  Amount  of  aid  (lbs)  transported  (per  trip)  from  depot  d  to  storage  facility  j  in 
scenario  s 


loop(d, 

loopG, 

loop(s, 

loop(b, 

Atd(d,j,s,b)  =  64000; 

); ); ); ); 


*  Amount  of  aid  (lbs)  transported  (per  trip)  from  storage  facility  j  to  disaster  area  i 
via  delivery  mode  k  in  scenario  s 

loop(i, 

loopG, 

loop(s, 

loop(b, 

Ats(i,j,'C17',s,b)  =  64000; 

Ats(i,j,'C130',s,b)  =  25600; 

Ats(i,j,’CH47’,s,b)  =  10000; 

); ); ); ); 

*  Inventory  holding  cost  (per  lb)  at  depot  d 
IHDcost(’Dl’)  =  unifonn(3,7); 

IHDcost('D2’)  =  unifonn(2,6); 

*  Inventory  holding  cost  (per  lb)  at  storage  facility  j 
IHScost('Fl’)  =  unifonn(l,  10); 

IHScost('F2')  =  unifonn(l,  10); 
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IHScost('F3’)  =  uniform(l,  10); 

IHScost('F4')  =  uniform(l,  10); 

IHScost('F5')  =  uniform(l,  10); 

IHScost('F6')  =  unifomi(l,  10); 

IHScost('F7')  =  uniform(l,  10); 

IHScost('F8’)  =  uniform(l,  10); 

IFIScost(T9’)  =  uniform)  1,  10); 

IHScost('FlO’)  =  unifonn(l,  10); 

*  Procurement  cost  (per  lb)  of  consumable  aid  from  vendor  v 
Pcost('Vl')  =  uniform(  10,20); 

Pcost('V2')  =  uniform(23,35); 

Pcost('V3’)  =  uniform(15,27); 

*  Inventory  holding  capacity  (lbs)  at  storage  facility  j 
Cap('Fl')  =  unifonn(2 1000000, 105000000); 

Cap('F2')  =  unifonn(2 1000000, 105000000); 

Cap('F3')  =  unifonn(2 1000000, 105000000); 

Cap('F4')  =  unifonn(2 1000000, 105000000); 

Cap('F5')  =  unifonn(2 1000000, 105000000); 

Cap('F6')  =  unifonn(2 1000000, 105000000); 

Cap('F7')  =  unifonn(2 1000000, 105000000); 

Cap('F8')  =  unifonn(2 1000000, 105000000); 

Cap('F9')  =  unifonn(2 1000000, 105000000); 

Cap('FlO’)  =  unifonn(2 1000000, 105000000); 

*  Inventory  holding  capacity  (lbs)  at  depot  d 
Cp(’Dl')  =  unifonn(250000000, 350000000); 

Cp('D2')  =  unifonn(250000000, 350000000); 

*  Fixed  cost  of  opening  storage  facility  j 
Fcost('Fl')  =  uniform)  100000,200000); 

Fcost('F2')  =  uniform)  100000,200000); 

Fcost('F3')  =  uniform)  100000,200000); 

Fcost('F4')  =  uniform)  100000,200000); 

Fcost('F5')  =  uniform)  100000,200000); 

Fcost('F6')  =  uniform)  100000,200000); 

Fcost('F7')  =  uniform)  100000,200000); 

Fcost('F8')  =  uniform)  100000,200000); 

Fcost('F9')  =  uniform)  100000,200000); 

Fcost('FlO')  =  uniform)  100000,200000); 

*  Maximum  number  of  trips  with  consumable  aid  from  depot  d  to  storage  facility 

j 

*  based  upon  a  set  replenish  time  target  of  120  hours  (5  days) 
loop(d, 
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loopCj, 

Ntd(d,j)  =  round(120  /  (2  *  r(d,j))); 

);); 

*  Maximum  number  of  trips  with  consumable  aid  from  storage  facility  j  to 
disaster  area  i 

*  via  aerial  delivery  mode  k;  based  upon  a  set  response  time  target  (eg.  72  hours) 
loop(i, 

loop(j, 

loop(k, 

Nts(i,j,k)  =  round(72  /  (2  *  t(i,j,k))); 

); ); ); 


*  Number  of  aerial  delivery  modes  k  pre-positioned  at  storage  facility  j 
Table  ADM(j,k) 

$ondelim 

$include  NumModes.csv 
$offdelim 


*  Establish  the  optimization  iterations 

Table  its(b,  s,  fact)  DOE  factors  for  each  scenario  and  iteration  combination 
$ondelim 

$include  Iterations. csv 
$offdelim 


* _ 

*  Decision  Variables  and  Goal  Deviation  Variables 

*  _ 

Binary  Variables 

X(j)  equals  1  if  "open"  storage  facility  j  or  0  otherwise 

Y(i,j)  equals  1  if  disaster  area  i  is  covered  by  storage  facility  j  or  0 

otherwise 


Integer  Variables 

Z(i,j,k)  number  of  trips  with  consumable  aid  from  storage  facility  j  to 
disaster  area  i  using  aerial  delivery  mode  k 

H(d,j)  number  of  trips  with  consumable  aid  from  depot  d  to  storage  facility 

j 


Positive  Variables 
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Q(d,v) 

Inv(j) 

In(d) 

pos(g,s,b) 

neg(g,s,b) 


amount  of  aid  (lbs)  to  purchase  from  vendor  v  for  storage  at  depot  d 
amount  of  aid  (lbs)  to  hold  in  inventory  at  storage  facility  j 
amount  of  aid  (lbs)  to  hold  in  inventory  at  depot  d 
positive  deviation  for  goal  g  in  scenario  s  and  iteration  b 
negative  deviation  for  goal  g  in  scenario  s  and  iteration  b 


Free  Variable 

M  objective  function  value  to  be  minimized 


Equations 

OBJ  define  the  overall  goal  programming  objective  function 

gc  1  goal  constraint  #  1 

gc2  goal  constraint  #2 

gc3  goal  constraint  #3 

gc4  goal  constraint  #4 

re  1  real  constraint  #  1 

rc2  real  constraint  #2 

rc3  real  constraint  #3 

rc4  real  constraint  #4 

rc5  real  constraint  #5 

rc6  real  constraint  #6 

rc7  real  constraint  #7 

rc8  real  constraint  #8 

rc9  real  constraint  #9 

re  1 0  real  constraint  #10 

re  1 1  real  constraint  #11 

re  1 2  real  constraint  #12 


* _ 

*  Objective  Function 

*  _ 

*  Minimize  the  sum  of  the  expected  weighted  goal  deviations  across  all  scenarios 
and  iterations 

OBJ..  M  =e=  sum((b,s),  its(b,  s,'prob')*(w(T)*pos(T,s,b)  +  w(’2')*pos(’2',s,b)  + 
w(’3')*neg(’3’,s,b)  +  w('4')*pos(’4’,s,b))); 

*  _ 

*  Goal  Constraints 

*  _ 

*  Minimize  the  total  response  time  (in  trip-hours)  for  each  scenario  and  iteration 
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gcl(s,b)..  sum((d,i,j,k),  t(i,j,k)*Z(i,j,k)  +  r(d,j)*H(d,j»  +  neg('l',s,b)  -  pos(’l’,s,b) 
=e=  0; 

*  Minimize  the  total  supply  chain  cost  for  each  scenario 
gc2(s,b)..  sum((d,i,j,k,v),  its(b,s,'SCDisrupf)*(TDcost(d,j)*H(d,j)  + 
TScost(i,j,k)*Z(i,j,k)  +  Fcost(j)*X(j)  +  IHDcost(d)*In(d)  +  IHScost(j)*Inv(j)  + 
Pcost(v)*Q(d,v)))  +  neg('2’,s,b)  -  pos(’2',s,b)  =e=  0; 

*  We  want  the  amount  of  consumable  aid  delivered  to  meet  the  amount  of 
expected  demand  for  each  scenario 

gc3(s,b)..  sum((i,j,k),  Ats(i  j ,k,s,b)*Z(i  j ,k)*AD M (j ,k))  +  neg(’3’,s,b)  -  pos(’3’,s,b) 
=e=  sum(i,  Dem(i,s,b))*its(b,s,'dishnp')*its(b,s,'DaysFood'); 

*  Minimize  the  amount  of  expected  demand  not  met  (shortage)  for  each  scenario 
gc4(s,b)..  neg(’3',s,b)  +  neg('4’,s,b)  -  pos('4',s,b)  =e=  0; 


* _ 

*  Real  Constraints 

*  _ 

*  Each  disaster  area  i  must  be  covered  by  one  storage  facility  j : 

rcl(i)..  sum(j,  Y(i,j))  =e=  1; 

*  If  disaster  area  i  is  covered  by  storage  facility  j,  then  storage  facility  j  must  be 
opened: 

rc2(j)..  sum(i,  Y(i,j))  =e=  X(j); 

*  The  amount  of  inventory  to  store  at  storage  facility  j  must  be  greater  than 

*  or  equal  to  the  total  product  of  the  number  of  trips,  the  amount  of 

*  consumable  aid  transported  per  trip,  and  the  number  of  assets: 
rc3(j,s,b)..  sum((i,k),  Ats(i,j,k,s,b)*Z(i,j,k)*ADM(j,k))  =1=  Inv(j); 

*  The  inventory  is  held  only  at  "opened"  storage  facilities,  and  the  amount  of 

*  inventory  kept  at  the  storage  facility  j  must  not  exceed  its  capacity: 

rc4(j)..  Inv(j)  =1=  Cap(j)*X(j); 

*  If  a  storage  facility  is  opened,  then  inventory  must  be  stored  there 
rc5(j)..  Inv(j)  =g=  X(j); 

*  The  amount  of  inventory  to  store  at  depot  d  must  be  greater  than  or  equal  to 

*  the  total  product  of  the  number  of  trips  and  the  amount  of  consumable  aid 

*  transported  per  trip: 

rc6(d,s,b)..  sum(j,  Atd(d,j,s,b)*H(d,j))  =1=  In(d); 

*  All  depots  are  assumed  to  be  open,  but  the  amount  of  inventory  kept  at  depot  d 
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*  must  not  exceed  its  capacity: 
rc7(d)..  In(d)  =1=  Cp(d); 


*  If  consumable  aid  is  purchased  from  the  vendors  for  storage  at  depot  d,  then 

*  there  is  inventory  for  storage  at  depot  d 
rc8(d)..  sum(v,  Q(d,v))  =1=  In(d); 

*  The  number  of  trips  with  consumable  aid  from  each  depot  d  to  the  each 

*  storage  facility  j  must  not  exceed  its  maximum: 
rc9(d,j)..  H(d,j)  =1=  Ntd(d,j)*X(j); 

*  The  number  of  trips  with  consumable  aid  from  each  storage  facility  j  to  each 

*  disaster  area  i  via  each  delivery  mode  k  must  not  exceed  its  maximum: 
rclO(i,j,k)..  Z(i,j,k)  =1=  Nts(i,j,k)*Y(i,j); 

*  The  total  amount  of  consumable  aid  shipped  from  the  depots  to  each  storage 
facility  j 

*  must  appropriately  backfill  the  total  amount  of  aid  delivered  to  the  disaster 
areas: 

rcl  l(j,s,b)..  sum(d,  Atd(d,j,s,b)*H(d,j))  =g=  sum((i,k), 
Ats(i,j,k,s,b)*Z(i,j,k)*ADM(j,k)); 

*  The  amount  procured  from  each  of  the  vendors  for  storage  at  all  of  the  depots 
must  be  greater 

*  than  or  equal  to  the  total  amount  of  aid  shipped  out  of  depot  d 
rcl2(v,s,b)..  sum(d,  Q(d,v))  =g=  sum((d,j),  Atd(d,j,s,b)*H(d,j)); 

*  We  set  the  bounds  on  the  decision  variables 

X. loG)=0; 

x.upO')=i; 

Y. lo(i,j)=0; 

Y. up(i,j)=l; 

Z. lo(i,j,k)=0; 

H.lo(d,j)=0; 

lnv.lo(j)=0; 

ln.lo(d)=0; 

Q.lo(d,v)=0; 

pos.lo(g,s,b)=0; 

neg.lo(g,s,b)=0; 

*  Setup  the  optimization  model 
model  mod  /all/; 

*  Suppress  the  number  of  rows  listed  to  0 
Option  limrow=0; 
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*  Suppress  the  number  of  columns  listed  to  0 
Option  limcol=0; 

*Sets  relative  optimality  tolerance  -  i.e.  no  tolerance  because  0  =  OPT  value 
Option  optcr=0.01; 

*  Solve  the  MILP  model  using  CPLEX  solver 
mod.OptFile=l; 

Option  MIP=cplex; 

SOLVE  mod  using  MIP  minimizing  M; 

*  Post-hoc  statistical  analysis 

display  X.l,  Y.l,  Z.l,  H.l,  Q.l,  Inv.l,  In.l,  pos.l,  neg.l; 

$ontcxt 

*===  Export  to  Excel  using  GDX  utilities 

*===  First  unload  to  GDX  fde  (occurs  during  execution  phase) 

execute  unload  " results l.gdx"  X.l  Y.l  Z.l  H.l  Q.l  Inv.l  In.l  pos.l  neg.l  t  r  TDcost 

TScost 


*===  Now  write  to  variable  levels  to  Excel  file  from  GDX 
*===  Display  Variables  Values  in  Excel 

execute  'gdxxrw.exe  results  l.gdx  o=resultsl.xls  var=X.l  rng=PrimalX!' 
execute  ’gdxxrw.exe  results  l.gdx  o=resultsl.xls  var=Y.l  rng=PrimalY!' 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=Z.l  mg=PrimalZ!’ 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=H.l  rng=PrimalH!' 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=Q.l  rng=PrimalQ!' 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=Inv.l  mg=PrimalInv!' 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=In.l  mg=PrimalIn!’ 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=pos.l  mg=Primalpos!' 
execute  'gdxxrw.exe  resultsl.gdx  o=resultsl.xls  var=neg.l  rng=Primalneg!' 


*===  Display  Parameter  Values  in  Excel 
execute  'gdxxrw.exe  resultsl.gdx  o=results  1  .xls 
execute  'gdxxrw.exe  resultsl.gdx  o=results  1  .xls 
execute  'gdxxrw.exe  resultsl.gdx  o=results  1  .xls 
execute  'gdxxrw.exe  resultsl.gdx  o=results  1  .xls 
$offtext 


par=t  rng=TravelTime_ijk!' 
par=r  rng=TravelTime_dj !’ 
par=TDcost  mg=TransCost_dj!' 
par=TScost  rng=TransCost_ijk!’ 


*===  Export  to  CSV  files 

file  facility  /FacilityOpen.csv/; 

facility. pc=5; 

put  facility; 

loop(O'), 

put  j.tl,  X.l(j)/; 

); 
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putclose; 

file  coverage  /DisasterCoverage.csv/; 

coverage.pc=5; 

put  coverage; 

loop((ij)$Y.l(ij), 

put  i.tl,  j.tl,  Y.l(i,j)/; 

); 

putclose; 

file  tripsl  /NumTripsDisaster.csv/; 

trips  l.pc=5; 

put  tripsl; 

loop((i,j,k), 

put  j.tl,  k.tl,  Z.l(i,j,k)/; 

); 

putclose; 

file  trips2  /NumTripsStaging.csv/; 

trips2.pc=5; 

put  trips2; 

l°°P((d,j), 

put  d.tl,  j.tl,  H.l(d,j)/; 

); 

putclose; 

file  procure  /ProcureVendor.csv/; 

procure.pc=5; 

put  procure; 

loop((d,v), 

put  d.tl,  v.tl,  Q.l(d,v)/; 

); 

putclose; 

file  invl  /InvStorage.csv/; 

invl.pc=5; 

put  invl; 

loop(O'), 

put  j.tl,  Inv.l(j)/; 

); 

putclose; 

file  inv2  /InvDepot.csv/; 

inv2.pc=5; 

put  inv2; 

loop((d), 


61 


put  d.tl,  In.  1(d)/; 

); 

putclose; 

fde  under  /UnderAchieve.csv/; 
under  .pc=5; 
put  under; 
loop((g,s,b), 

put  g.tl,  s.tl,  b.tl,  neg.l(g,s,b)/; 

); 

putclose; 

fde  over  /OverAchieve.csv/; 

over.pc=5; 

put  over; 

loop((g,s,b), 

put  g.tl,  s.tl,  b.tl,  pos.l(g,s,b)/; 

); 

putclose; 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ANG 

Air  National  Guard 

Cl 

confidence  interval 

DOD 

Department  of  Defense 

DOE 

Design  of  Experiments 

DOS 

Department  of  State 

DRAP 

Disaster  Relief  Airlift  Planner 

EM-DAT 

Emergency  Events  Database 

GAMS 

General  Algebraic  Modeling  System 

GHz 

gigahertz 

HA/DR 

Humanitarian  Assistance/Disaster  Relief 

LB 

lower  bound 

MCDA 

Multiple  Criteria  Decision  Analysis 

MCLB 

Marine  Corps  Logistics  Base 

NSRDEC 

US  Army  Natick  Soldier  Research,  Development  and 
Engineering  Center 

RDECOM 

US  Army  Research,  Development  and  Engineering 

Command 

SAA 

Sample  Average  Approximation 

UB 

upper  bound 

USAID 

US  Agency  for  International  Development 

USTRANSCOM 

US  Transportation  Command 

VBA 

Visual  Basic  for  Applications 

VTD 

Vehicle  Technology  Directorate 

WGSS 

Within-Group  Sum  of  Squares 
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